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SU(2)  Hamiltonian  lattice  gauge  theory  is  used  as  a  theoretical  model  to  test 
new  computational  methods  developed  in  this  dissertation  and  to  study  the 
possible  features  exhibited  by  Hamiltonian  lattice  QCD. 

The  basic  method  used  here  to  study  SU(2)  Hamiltonian  lattice  gauge 
theory  is  a  variational  calculation,  which  enables  non-perturbative  calculations 
to  be  made  from  the  strong-coupling  region  down  to  the  scaling  region. 
Expectation  values  of  physical  observables  are  calculated  using  a  completely 
analytical  method  requiring  approximation  of  the  localization  of  the  magnetic 
fluctuations.  By  examining  the  analytical  calculation  method  in  detail,  the 
nature  and  general  rules  for  applications  of  the  method  are  established. 

The  accuracy  of  this  method  as  well  as  the  potential  speed  and  accuracy  of 
a  new  computer  algorithm,  the  Biased-Selection  Monte  Carlo  method,  are 
demonstrated  by  reproducing  numerically  many  of  the  values  of  certain  physical 
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observables  without  making  the  approximation  required  by  the  analytical 
method.  The  results  demonstrate  that  this  new  algorithm  is  accurate  and  also  is 
reasonable  in  its  consumption  of  computer  time. 

The  physical  observables  studied  are  the  static  quark  potential  and  the 
glueball  mass  of  SU(2)  Hamiltonian  lattice  gauge  theory.  The  product  of  a  simple 
plaquette  wavefunction  is  used  as  the  variational  ground  state.  From  the  heavy 
quark  potential,  the  rough  transition  point  is  located  at  ^c  =  l/2gc  =  0.26.  The 

string  tension  also  can  be  extracted  from  the  heavy  quark  potential,  the  string 

tension  exhibits  scaling  behavior  which  is  consistent  with  the  expected  asymptotic 

scaling  line.  The  scaling  behavior  of  the  string  tension  implies  that 
=  0.009  Jo ,  o  being  the  string  tension. 

The  lowest  excited  state,  including  correlations  of  up  to  two  adjacent 
plaquettes,  is  found  by  using  a  projection  operator  to  remove  the  ground  state. 
This  enables  calculation  of  the  glueball  mass.  The  result  exhibits  an  obvious 
tendency  toward  seating,  but  the  glueball  mass  in  this  approximation  does  not 
scale  exacdy  along  the  expected  asymptotic  scaling  line.  However,  a  scaling 
window  is  observed  in  the  ratio  of  glueball  mass  to  the  square  root  of  suing 
tension,  and  it  yields 

M„  =  2.86  Jo 
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CHAPTER  1 
INTRODUCTION  AND  OVERVIEW 


Quantum  chromodynamics  (QCD)  now  is  accepted  widely  as  an  underlying 
theory  of  strong  interaction.  This  theory  has  two  quite  different  aspects.  On 
one  hand,  because  of  the  property  of  asymptotic  freedomO)  possessed  by  the 
non-abelian  gauge  theories,  at  short  distance  the  effective  coupling  constant  of 
the  quark-gluon  interaction  becomes  small.  The  interaction  can  then  be  treated 
perturbatively  and  the  theory  becomes  relatively  simple.  On  the  other  hand, 
quarks  within  hadrons  are  bounded  and  inseparable.  In  this  situation,  the 
dynamical  mechanism  becomes  very  complicated.  The  complexity  can  be  seen 
most  easily  from  the  fact  that  the  proton  is  built  out  of  three  quarks  (u,  u  and  d) 
which  are  much  tighter  than  the  proton.  In  order  to  demonstrate  that  QCD  is  a 
comprehensive  theory,  the  complexity  must  be  overcome  and  large-scale 
dynamics  must  be  included. 

To  this  end,  lattice  gauge  theory  was  introduced  by  Wilson(2)  as  a 
theoretical  aid  and  computation  tool  to  study  the  large-scale  properties  of  QCD. 
Once  a  field  theory  has  been  formulated  on  a  lattice,  the  problem  becomes  one 
of  statistical  mechanica,  i.e.,  reduces  in  a  sense  to  the  evolution  of  the 
multidimensional  integral  known  as  the  partion  function.  In  this  work,  we  used 
both  analytical  and  Biased-Selection  Monte  Carlo  methods  to  evaluate  this 
integral  at  various  values  of  the  coupling  constant  accurately  enough  to  extract 
information  about  the  physical  contents  of  the  system.  To  understand  the 
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underlying  physics,  we  have  to  probe  the  different  phases  of  this  system,  and 
understand  the  nature  of  its  critical  points.  In  particular,  if  the  transition 
between  two  different  phases  of  the  system  is  continuous,  we  can  obtain  a  well 
defined  continuum  limit  by  approaching  the  critical  point  of  the  theory.  The 
phases  of  the  statistical  mechanics  system  can  be  obtained  in  a  variety  of 
methods,  such  as  duality  transformation,  high  or  low  temperature  expansion, 
renormalization  group  technique,  etc.  For  example,  the  high  temperature 
expansion  yields  the  following  form  of  the  two  point  correlation  function  of  3- 
dimension  Ising  model  (3) 

<  Sk  Sn  >  ~  exp(  -  mr ) 

as  r  =  Ik  -  nl->  °°i  here  m  is  a  constant  Using  duality  transformation,  one  can 
show  that  in  the  low  temperature  regime  the  correlation  function  above  behaves 
as 

t-  n  a  /  area(c) . 
rc  =  <gV  ~ exP  (  f^> 

where      are  the  dual  gauge  variables,  area(c)  is  the  minimal  area  enclosed  by 

c(a  closed  curve  over  the  gauge  variables  on  the  lattice  links),  and  T  is  the 
temperature.  Analogously,  Wilson  introduced  the  Wilson  loop(2)  as  an  order 
parameter  for  the  nonabelian  gauge  theories.  It  is  defined  as 
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<  W(c)  >  =  Tr P exp  ( i |      dx* )  > 

c 

=  JdA  e"S(A)  Tr  P  exp  ( i  |  A^  dx^ ). 

c  (1-1) 


Using  perturbative  techniques  in  the  weak  coupling  region,  one  can  show  that 
these  nonabelian  gauge  theories  are  asymptotically  free.  In  addition,  if  the 
above  Wilson  loop  obeys  an  area  law 


<  W(c)  >  ~  exp(-area(c)) 


(1-2) 


one  can  show  that  the  theory  linearly  confines.  Since  the  introduction  of  lattice 
gauge  theory,  many  qualitative  analyses  and  numerical  simulations  have  been 
performed,  all  of  which(^)  indicate  that  the  color  confinement  coexists  with 
asymptotic  freedom  in  QCD,  and  there  is  no  phase  transition  between  them. 

Because  the  formulation  of  Euclidean  lattice  gauge  theory  preserves 
reflection  positivity  and  admits  a  self-adjoint  transfer  matrix  (see  Chapter  2), 
there  is  also  an  alternative  formalism,  the  Kogut-Susskind  Hamiltonian(5)  of 
SU(N)  lattice  gauge  theory,  which  is  useful  in  studying  the  low-lying  spectrum 
of  the  theory.  In  this  dissertation,  all  of  the  analytical  calculations  and  the 
numerical  simulations  are  done  in  the  Hamiltonian  formalism.  There  is  a  basic 
philosophy  behind  working  in  the  Hamiltonian  formulation  to  deal  with  a 
complicated  theory  such  as  QCD.  For  example,  consider  the  relation  between 
the  well-analyzed  theory  of  Quantum  Electrodynamics  (QED)  and  the 
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extensively  studied  phenomena  of  superconductivity.  In  principle,  it  is  possible 
to  start  from  the  first  principle,  QED,  and  the  existence  of  nuclei  and  electrons, 
first  to  establish  the  crystal  structure  of  solids,  then  from  that  the  electron- 
phonon  interaction,  then  the  BCS  theory ,(6)  and  finally  superconductivity.  In 
practice,  however,  no  one  has  been  able  to  succeed  even  in  the  first  step.  One 
need  not  be  overly  surprised  if  one  finally  has  to  work  in  a  phenomenological 
model  of  hadron  physics  because  of  the  complexity  in  QCD.  The  crucial  point 
of  working  in  the  Hamiltonian  formalism  is  to  determine  the  important  field 
configurations  which  should  be  included  in  the  wavef unction.  The  information 
thus  obtained  will  give  a  clear  physical  picture  to  the  large-scale  dynamical 
mechanism  of  QCD  and  will  provide  the  necessary  ingredients  for  constructing 
a  phenomenological  model,  which  will  go  beyond  the  rather  simplified  existing 
models  such  as  the  MITC7)  and  the  SLAC(8)  bag  models. 

To  perform  calculations  in  lattice  gauge  theory,  there  are  two  approaches 
currently  used  by  most  of  the  theorists  in  the  field: 

(a)    Analytical  calculations  by  the  use  of  strong  coupling  expansion:  In 
the  limit  of  large  lattice  spacing,  the  theory  is  in  a  strong  coupling 
region,exhibits  confinement,  and  can  be  calculated  exactly.  If  there  are  no 
phase  changes,  perturbation  theory  and  the  techniques  of  Pade 
approximation^)  allow  calculations  also  in  the  weak  coupling  region.  In  this 
way,  small  lattice  spacing  and  continuum  physics  can  be  reached.  This  method 
is  widely  used  both  in  Lagrangian  and  Hamiltonian  formalisms.  However,  the 
success  of  this  approach  depends  critically  on  the  absence  of  any  phase 
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transition  (including  roughening  transition),  which  will  prevent  the  continuation 
to  the  weak-coupling  regime  where  the  continuum  field  theory  is  recovered. 

(b)    Numerical  simulations  by  the  use  of  Monte  Carlo  method:  With  this 
approach,  the  Monte  Carlo  method  is  applied  to  the  lattice  gauge  system  to 
evaluate  expectation  values.  The  idea  of  the  method  is  to  replace  the 
computation  of  the  expectation  value  of  a  quantity  9 


<9>  = 


(1-3) 


by  an  average  over  gauge  configurations  {U}i 

Xe[{u)i] 


(1-4) 


where  the  configurations {U}i  are  distributed  according  to  the  Boltzmann  weight 
e~S.  In  the  Hamiltonian  formalism,  e*S[Uk]  js  replaced  by  y2  [Uk]  where  y 
[Uk]  is  the  wavefunction.  To  replace  eq.  (1-3)  by  eq.  (1-4),  it  is  necessary  to 
single  out  the  important  configurations.  This  is  the  technique  called 
"importance  sampling. "(1°)  The  most  well-known  and  widely  used  one  is  the 
Markov  chain/10)  The  simple  Markov  chain  has  large  computer  time 
requirement,  and  bad  error  estimates.  It  is  therefore  of  interest  to  investigate  a 
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computer  algorithm  which  is  potentially  more  economical  in  the  consumption  of 
computer  time  and  also  gives  accurate  error  estimates. 

The  basic  computational  method  used  in  this  dissertation  is  the  variational 
method.  Starting  with  a  trial  wavefunction  based  on  physical  properties  of  the 
system,  calculations  are  carried  out  from  the  strong  coupling  region  down  to  the 
scaling  region  by  varying  variational  parameters  as  a  function  of  coupling 
constant  Finally,  the  result  in  the  scaling  region  is  compared  with  the  expected 
asymptotic  scaling  line  to  determine  if  the  continuum  physics  is  recovered.  To 
carry  out  this  procedure  analytically,  in  Chapter  3.2,  a  systematic  way  of 
performing  the  calculations  motivated  by  Batrouni's  mean  plaquette  method^  *) 
and  Robson's  plaquette- space  integration^)  js  developed-  In  Chapter  6.2  a 
new  computer  algorithm,  called  the  Biased-Selection  Monte  Carlo  method,  is 
introduced.  This  algorithm  was  first  applied  to  a  hard  disk  system  of  statistical 
mechanics  by  Coldwell.03)  As  a  demonstration,  the  Biased  Selection  Monte 
Carlo  method  is  used  to  estimate  some  physical  quantities.  The  results  indicate 
that  this  new  algorithm  can  overcome  problems  met  by  other  available 
techniques. 

Following  the  above  strategy,  using  the  analytical  method  developed  in 
Chapter  3  two  properties  of  SU(2)  Hamiltonian  lattice  gauge  theory  (gluon 
sector  only)  are  studied  in  detail  in  the  dissertation:  The  field  configuration  of 
the  ground  state  and  the  first  excited  state.  For  the  ground  state,  the  ground 
state  energy  per  plaquette  and  static  quark  potential  is  studied  in  Chapter  4;  for 
the  first  excited  state,  the  glueball  mass  is  studied  in  Chapter  5.  The  gauge 
group  SU(2)  was  studied  because  of  the  simplicity  of  its  group  structure  and  yet 


similar  basic  physical  contents  of  SU(3).  But,  in  some  ways,  the  case  of  SU(2) 
appears  to  be  more  complicated  than  that  of  the  more  realistic  gauge  group 
SU(3)  because  of  correlation  effects.  The  results  obtained  here  are  helpful  in 
establishing  rigorous  computational  methods,  and  in  understanding  the 
underlying  physics. 


CHAPTER  2 
FORMULATION 


In  1974,  Wilson  (2)  introduced  an  elegant  formalism  of  gauge  theories  on  a 
space-time  lattice.  The  Minkowskian  space-time  first  is  transformed  through  a 
Wick  rotation  into  an  Euclidian  space-time.  The  Euclidian  space-time  then  is 
discretized  to  act  as  an  ultraviolet  cut-off,  and  the  theory,  which  is  formulated 
on  a  discretized  Euclidean  space-time  lattice,  becomes  equivalent  to  a  statistical 
mechanics  problem.  The  most  remarkable  property  of  this  formalism  is  that 
local  gauge  invariance  always  remains  exact.  Besides  Wilson's  formalism, 
there  is  an  alternative  one,  first  advocated  by  Kogut  and  Susskind.(5)  In  the 
Kogut-Susskind  Hamiltonian  formalism,(5)  the  time  variable  is  kept  continuous 
and  the  space  variable  is  discretized.  Working  in  the  temporal  gauge  Aq  =  0, 
they  define  a  Hamiltonian  which  is  a  function  of  the  space  components  of  the 
gauge  fields  and  a  set  of  the  conjugate  momenta.  Instead  of  being  analogous  to 
statistical  mechanics  in  Wilson's  formalism,  a  conventional  quantum  mechanical 
system  with  a  well-defined  Hamiltonian  is  dealt  with.  In  Wilson's  formalism, 
strong  coupling  expansion  is  analogous  to  high  temperature  expansion  in 
statistical  mechanics.  The  Hamiltonian  formalism  also  permits  a  strong 
coupling  expansion,  which  is  now  an  application  of  quantum  mechanical 
perturbation  theory. 

In  this  chapter,  the  Wilson  action  will  be  constructed  first;  then  the  Kogut- 
Susskind  Hamiltonian  will  be  derived  from  the  Wilson  action.  Some 
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discussions  and  comments  are  made  for  both  Lagrangian  and  Hamiltonian 
formalisms.  Finally,  the  scaling  behavior  of  lattice  gauge  theory  and  the 
continuum  physics  will  be  discussed.  The  information  contained  in  this  chapter 
is  not  completely  new  for  the  researchers  in  this  field.  However,  here  the 
author  has  attempted  to  clear  up  some  clouds  ,which  once  puzzled  the  author, 
about  the  fundamental  formulations  of  lattice  gauge  theories. 

2.1.  Lagrangian  Lattice  Gauge  Theory 
The  action  is  a  mathematical  construction  with  two  important  properties: 

(a)  According  to  the  Hamiltonian  principle,  the  classical  equations  of 
motion  can  be  obtained  by  extremizing  the  action. 

(b)  According  to  the  Noether  theorem,  the  symmetries  of  the  action 
correspond  to  the  various  conserved  quantities  in  the  physical  system. 

Therefore,  a  guide  for  constructing  an  action  is  usually  the  symmetry 
respected  by  the  physical  system.  Today,  physicists  believe  that  the  symmetry 
associated  with  the  Poincare  group,  i.e ,  the  inhomogeneous  Lorentz  group,  is 
respected  by  all  of  the  physical  systems.  In  addition,  the  spectacular  results 
obtained  in  QED,  Quantum  Flavordynamics  (QFD)  and  QCD  make  gauge 
symmetry  indispensable  for  constructing  an  action.  However,  even  if  the 
symmetries  of  the  physical  system  are  specified  there  still  are  unlimited  ways  to 
construct  an  action.  In  continuous  theory,  two  properties  usually  are  used  to 
single  out  a  particular  action: 

(a)    Renormalizability,  which  makes  the  ultraviolet  and  the  infrared 
infinities  manageable. 
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(b)  Simplicity,  as  long  as  the  dynamics  represented  by  the  action  are  non- 
trivial. 

The  action  constructed  in  this  way  is  called  the  "minimal  action."  In  lattice 
formulation  of  field  theory,  the  space-time  symmetry  of  the  Poincare  group  is 
discretized,  but  an  effort  is  made  here  to  maintain  the  exact  gauge  symmetry. 
Because  of  the  cut-off  introduced  by  the  space-time  lattice,  ultraviolet  infinities 
can  be  ignored,  but  the  action  is  constructed  to  be  as  simple  as  possible.  This  is 
exactly  how  the  Wilson  action  is  constructed. 
2.1.1.  The  Wilson  Action  of  SU(N)  Lattice  Gauge  Theory 

Considering  a  four-dimensional  Euclidean  space,  a  hypercubical  lattice 
with  spacing,  a,  is  introduced.  Links  are  defined  as  the  bonds  connecting  each 
pair  of  the  nearest  neighbor  sites.  There  are  eight  links  emerging  from  each 

site.  An  elementary  square,  which  is  enclosed  by  four  links,  is  called  a 
plaquette.  A  link  pointed  from  site  n  along  the  direction  fl  to  the  nearest 
neighbor  site  is  associated  with  an  independent  group  element  Uu.(n) 

iyn)E  SU(N)  (2.1} 

On  transversing  a  link  in  the  opposite  direction,  one  should  obtain  the  inversed 
element 


U.jl(n  +  A)  =  Ujl-1(n) 


(2-2) 
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The  U. .  (n^  matrices  are  defined  as 


(n)  =  exp(i  gD  aA ^  (n  +  ^  (1))  (2-3) 


where  g0  is  the  bare  coupling  constant.  The  vector  potential  A|j.  is  the  matrix 


AM-  "  ?  K  T&  (2-4) 


where  Ta  are  generators  of  the  Lie  algebra  of  gauge  group  SU(N).  The 
coordinate  associated  with  Ajj,  should  be  in  the  vicinity  of  the  link  in  question. 
For  convenience,  the  coordinate  is  assumed  to  be  in  the  middle  of  a  given  link. 

To  preserve  local  gauge  invariance,  the  action  should  consist  of  various 
closed  loops.  For  simplicity,  Wilson(2)  proposed  that  the  action  should  be  a 
sum  over  all  plaquettes  of  the  lattice 

S=rTSP 

P  (2-5) 
where  the  subscript  P  denotes  plaquette.  The  action  of  each  plaquette  Sp  is  the 
trace  of  the  product  of  the  group  elements  associated  with  the  four  links  which 
form  the  surface  boundaries  of  the  plaquette,  as  shown  in  Figure  1. 
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n  +  0 


U+(n+0) 


U+  (n) 


n 


Uv(n  +  &) 


n  +  ft 


Figure  1 .  A  plaquette  Up  with  Up  =      (n)  Uy  (n  +  (I)  U  *  (n  +  0)  Uy  +  (n) 


It  also  is  required  here  that  when  the  group  element  of  a  plaquette  approaches 
identity,  the  action  of  a  plaquette  should  approach  zero.  Therefore,  the  action  of 
a  plaquette  for  the  gauge  group  SU(N)  can  be  written  as 


Sp  =  P[l-I-RetrUp] 


(2-6) 


and,  for  a  plaquette  laying  in  the  ( \i,\> )  plane, 
Up  =  iyn)Uv(n+A)U*  (n+0)Uv  +  (n) 


(2-7) 
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The  normalization  factor  (3  will  be  defined  later.  The  trace  in  eq.  (2-6)  can  be 

taken  in  any  representation,  but  for  simplicity  only  the  fundamental 
representation  is  considered. 

To  demonstrate  that  the  action  of  eq.  (2-5)  reduces  to  the  usual  Yang-Mills 
theory  in  the  continuum  limit,  the  plaquette-action,  eq.  (2-6),  is  expanded  with 
respect  to  the  center  position  of  the  plaquette  in  terms  of  the  powers  of  lattice 
spacing,  a.  Consider,  for  example,  a  plaquette  lying  in  the  (1,2)  plane  with  the 
center  at  n^  and  with  counterclockwise  orientation.  The  corresponding  action  is 


Sp  =  P{l-iRetr[exp(ig0aA1(n 


)] 


•  exp(-igoaA20y-i-a8  ))]} 


(2-8) 


Taylor  expanding  the  vector  potential  about  the  center  of  the  plaquette  gives,  up 
to  the  order  0(a2): 
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Sp  =  P  { 1  - 1-  Re  tr  [exp(  i  gQ  aA,  (ty)  - 1  gQ  a2  32  A,  (iy) 

•  exp  ( i  gQ  aA2  (iy)  +  y  gQ  a^  A2  (n  ^  ) 

•  exp  ( -i  gQ  aAj  (ry)  -  ±  gQ  a2  B2  Aj  ) 

•  exp  ( -i  gQ  aA2       +  ^  gQ  a2  3j  A2  (ly) )] } 

Use  of  the  Baker-Hausdorff  formula 


exp  (x)  exp  (y)  =  exp  (x  +  y  +  i-  [x,y]  +  ...) 

2  (2-10) 


(2-9) 


yields 


Sp  =  p  [1  -  JL  Re  Tr  exp  (i  g0  a2  Qfo  -  92  Aj  +  ig0  [A^l  ) )  ] 
=  P[l-iReTrexp(ig0a2F12)] 


P  „2_4^_^  ,2 

(2-11) 


2N  goa 


Sumrning  over  all  of  the  plaquettes  gives 
=  JL J  IV    1^  ,2 


S  =  2^go  Jd44Tr(F^)2  (2"12) 
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with  the  factor  of  one-half  inside  the  integral  coming  from  the  symmetry  under 
\ix>  interchange.  If  p  is  defined  as 


P=  2N 


go  (2-13) 


the  Wilson  action  reduces  to  the  usual  Yang-Mills  theory  in  the  limit  a  — >  0. 
2.1.2.  Discussions  and  Comments 

From  the  view  of  constructive  quantum  field  theory,  Wilson's  lattice 
formalism  of  gauge  theories  has  the  following  properties: 

(a)  The  lattice  regularization  preserves  local  gauge  invariance:  When 
space-time  is  discretized,  a  cut-off  is  introduced  into  the  theory  naturally.  The 
introduced  cut-off  is  proportional  to  1/a  in  momentum  space.  Usually  a  cut-off 
will  break  local  gauge  invariance,  as  in  the  well-known  case  of  the  Pauli-Villars 
regularization.(14)  But  as  part  of  the  theory,  local  gauge  invariance  is 
necessary  in  establishing  both  ultraviolet  and  infrared  properties  of  the  theory. 
This  is  exactly  the  reason  that  use  of  the  dimensional-regularizationC^)  in  the 
continuum  gauge  theories  is  preferred. 

(b)  The  formulation  preserves  reflection  positivity  and  admits  a 
self-adjoint  transfer  matrix:^)  Because  of  reflection  positivity,  a  positive 
definite  Hilbert  space  of  physical  states  can  be  reconstructed;  and  because  of  the 
existence  of  a  self-adjoint  transfer  matrix,  a  self-adjoint  Hamiltonian  also  can  be 
defined.  These  two  characteristics  guarantee  that  the  formulation  specifies  an 
Euclidean  quantum  field  theory. 
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(c)  The  formulation  reduces  to  the  usual  Yang-Mills  theory  in  the  limit  of 
lattice  spacing  approaching  zero. 

Along  with  these  remarkable  properties,  there  are  some  problems: 

(a)  The  formulation  severely  mutilates  the  space-time  symmetry  of  the 
Poincare  group:  This  problem  reflects  the  fact  that  troubles  are  encountered  in 
trying  to  consistently  include  fermions  in  the  formulation.(^)  For  the  gauge 
fields,  the  Lorentz  properties  appear  to  be  irrelevant  to  the  lattice  formulation. 

(b)  The  formulation  obscures  the  topological  properties  of  gauge  fields: 
On  a  lattice,  continuity  in  space-time  is  lost  and  any  lattice  gauge  field  U|j.(n) 
can  be  deformed  continuously  to  the  trivial  field  U(i(n)  =  1,  where  apparently 
there  is  no  topological  structure. 

As  far  as  these  problems  are  concerned,  it  is  noted  that  in  the  weak- 
coupling  region  the  space-time  symmetry  will  be  restored  eventually  and  gauge 
fields  again  will  fall  into  disconnected  topological  classes.  Therefore, 
topological  interpretations  can  be  given  to  some  field  configurations  according 
to  the  topological  classes  with  which  they  associate  in  the  weak-coupling 
region.  Some  work  has  been  done  in  this  direction.^ 8) 

Finally,  it  is  desirable  to  comment  on  the  Wilson  action.  In  Wilson's 
action,  usually  the  fundamental  representation  of  a  gauge  group  is  used.  This  is 
convenient,  but  there  is  no  natural  reason  for  making  this  choice.  On  the  other 
hand,  in  the  previous  section  it  was  showed  that  Wilson's  action  reduces  to  the 
standard  Yang-Mills  action  in  the  naive  continuum  limit.  The  latter  property  is  a 
necessary  requirement  for  any  lattice  regularization  of  the  continuum  action. 
However,  it  is  not  sufficient  to  uniquely  determine  the  form  of  the  action. 
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Therefore  it  is  possible  to  have  various  actions,  all  with  the  same  continuum 
limit.  The  question  becomes  how  they  approach  the  continuum  limit  or,  if  the 
action  approaches  the  continuum  limit,  in  a  way  smooth  enough  to  avoid  any 
artificial  phase  transition.  Based  on  this  consideration,  Manton(19)  proposed 
another  action  which  shares  many  of  the  properties  of  the  Wilson  action.  In 
Manton's  approach,  the  action  is  constructed  in  terms  of  the  shortest  distance  in 
the  group  space 

SM  (U)  =  J?  I(d(U(p),l))2 

8    P  (2-14) 

where  d(U(p),  1)  denotes  the  geodesic  from  U  (p)  to  the  unit  matrix.  The 
action  here  is  independent  of  the  group  representation,  and  approaches  the 
continuum  limits  more  smoothly  than  the  Wilson  action  does/1^)  Lang  et 
alX2^)  used  Manton's  action  to  perform  calculations  for  U(N)  and  SU(N) 
lattice  gauge  theories  in  two  dimensions.  Their  results  showed  that  there  is  no 
sign  of  phase  transition  at  all,  even  for  the  U(oo)  theory  in  which  Gross  and 
Witten(21)  observed  a  third-order  phase  transition  by  using  Wilson's  action. 
However,  as  stressed  by  Kogut(22)  and  others,  the  detail  form  of  the  lattice 
action  should  become  irrelevant  in  the  continuum  limit  While  most  of  the 
researchers  in  this  field,  including  this  author,  still  use  the  Wilson  action,  it  is 
necessary  to  be  cautious  about  any  transition  exhibited  by  the  theory.  Monte 
Carlo  studies  of  the  Wilson  action  show  a  rapid  but  smooth  crossover  from 
ordered  to  random  behavior  for  SU(2)  and  SU(3),(4)  and  a  first-order  phase 
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transition  for  SU(N)  with  N  >  4  .(23)  xhe  \asi  generally  is  believed  to  be  an 
artifact  of  the  Wilson  action,  rather  than  a  transition  to  a  new  phase  of  the 
continuum  theory. 

(24) 

Therefore,  in  regard  to  SU(2)  and  SU(3),  it  is  safe  to 

use  the  Wilson  action. 

2.2.  Hamiltonian  Lattice  Gauge  Theory 

The  Kogut-Susskind  Hamiltonian  (5)  of  SU(N)  lattice  gauge  theory  was 
constructed  originally  according  to  the  gauge  symmetry.  In  this  section  for  the 
first  time  the  Hamiltonian  is  rigorously  derived  directly  from  the  Wilson  action 
without  using  the  transfer  matrix,  (using  the  transfer  matrix  to  derive  the  Kogut- 
Susskind  Hamiltonian  rigorously  can  be  founded  in  ref.(25)),  so  that  the 
equivalence  can  be  seen  between  these  two  formalisms.  To  derive  the 
Hamiltonian  from  the  Wilson  action,  the  rules  of  canonical  quantization  are 
followed:  The  Wilson  action  first  is  cast  in  Hamiltonian  form,  and  then 
canonical  commutation  relations  between  the  coordinates  and  their 
corresponding  conjugate  momenta  are  imposed.  Also  discussed  are  the  typical 
problems  appearing  in  the  canonical  quantization  of  continuous  nonabelian 
gauge  theory(26)  and  the  corresponding  status  of  lattice  gauge  theory. 
2.2.1.  Koeut-Susskind  Hamiltonian  of  SU(N)  Lattice  Gauge  Theory 

The  Wilson  action,  up  to  a  constant  term,  is 


S=-VsftrUD  +  Ii.c.} 


(2-15) 
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Since  it  is  intended  to  consider  the  continuous  time  limit  of  the  Wilson  action, 
the  space-time  lattice  is  deformed  by  introducing  different  lattice  spacings,  a  and 
ao,  for  space  and  time.  The  coupling  constants  on  spacelike  plaquettes  and 
timelike  plaquettes  no  longer  are  equal  in  the  Wilson  action 


S=^Z  {trU   +  h.c.}  +  -i  I   {trUp+  h.c.} 

1     1  s      s  (2-16) 


where  Pt  denotes  the  timelike  plaquette,  Ps  is  the  spacelike  plaquette,  gs  is  the 


bare  coupling  constant  for  the  spacelike  plaquette,  and  gt  is  that  for  the  timelike 
plaquette.  In< 

eq.  (2-16)  by 


2  1 

plaquette.  In  order  to  obtain  a  proper  classical  limit,  gt  andg~  are  replaced  in 


1    -4  a 


g2       g?aQ  (2-17a) 


and 


1     .  ao 


g2       g?  a  (2-17b) 


and  the  Wilson  action  becomes 


a 

S=-^—  I  {trU  +  h.c.)  +-^—  S  (trU    +  h.c.}  . 

gt  aD  Pt       pt  gs2a    Ps       Ps  (2-18) 
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In  the  limit  of  small  lattice  spacing,  i.e.,  a  ->  0  and  ao  ->  0,  we  have 


9  4  2 
trU_  +  h.c.  =  -  g  L  a  tr(F.)  +  constant 

Ps  s  1J  (2- 19a) 


and 


trLL  +  h.c.  =  -  g?  a2  a2  tr(F  .)2  +  constant 

pt  1       0      01  (2-19b) 


Here  i,j  =  1,2,  or  3.  Substituting  eqs.  (2-19a)  and  (2-19b)  into  eq.  (2-18),  a  proper 
classical  limit  is  obtained: 


S  ->^-Jd3xdtTr(F)m)2 


LU) 

a->U~" 

a0-^0  (2-20) 


To  proceed  toward  the  Hamiltonian  formahsm,  we  now  go  to  the  temporal  gauge 
by  setting  all  timelike  links  equal  to  identity.  The  timelike  plaquette,  as  shown  in 
Figure  2,  can  be  written  as 
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U+(t+l) 


U=  l\f 


A  U=  1 


U(t) 


►  x,  y,  z 

Figure  2.  A  time-like  plaquette  in  the  temporal  gauge. 


trUp  +  h.c.  =  tr(U(t)U+(t  +  1)  +  h.c.) 


=  -  tr{(  U+(t+l)  -  U+(t) )( U(t+1)  -  U(t) )}  +  constant 


(2-21) 


Substituting  eq.  (2-21)  into  eq.  (2-18), 


S  =  aQ{--a-  Str( 


Uf(t  +  1)  -  Uf(t)  U(t+1)-U(t) 


a„ 


ag^Ps  ps 


~^a0  {4r  I  trtlttI  +  -L 


ao^° 


0  g~  links  a  g?  ~Ps 


Itr(Un  +h.c.)}. 


(2-22) 
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From  eq.  (2-22),  in  the  continuum  time  the  Lagrangian  of  the  system  can  be  defined 


where  the  summations  over  links  and  plaquettes  only  refer  to  those  in  three 
spatial  dimensions.  The  action  then  is 


Now  it  is  possible  to  construct  the  Hamiltonian  through  the  canonical  procedure 


Note  that  while  working  in  the  temporal  gauge,  the  Lagrangian  L(t)  and  the 

Hamiltonian  H  still  have  the  invariance  of  time-independent  local  gauge 
rotation.  To  quantize  the  system  of  eq.  (2-25),      is  eliminated  in  favor  of  the 
generators  of  time-independent  local  gauge  rotation.  If  an  mfinitesimal  local 
gauge  rotation  is  made  at  site  n,  the  resulting  transformation  is 


(2-23) 


S=  dtL(t) 


(2-24) 


(2-25) 
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Uk  (n)  ->      (n)  =  Uk  (n)  +  i  £a  (n)  £  Uk  (n)  (2  26a) 


However,  because  links  are  orientated  and 


ut(„)  =  Uk(n  +  £)  (M6b) 


it  also  is  possible  to  have  the  transformation 


u  k  („)  -*  u;k  w = u  k  <„)  -  i  xa  w  u  k  („)  f 


a  ~a 

The  corresponding  two  sets  of  generators,  E^")      E^C0)'  satisfy 


[E^(n),U  (m)]=i^U  (n)  8  8 

k       J  2    k      j,k  n,m  (2-27a) 


t^(n),U.(m)]=-Uk(nD|  6  ^ 


[^(n),Ej,(m)]=ifabc  E^(n)8.k8nm 


(2-27b) 


(2-27c) 


and 


[E^n),^)]  =ifabc^(n)5jJc5nm 


(2-27d) 
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where  f300  are  structure  constants  of  the  gauge  group  and 


rl!  i1=ifabcx^ 

1  2  '  2  J  2    *  (2-28) 


According  to  the  Noether  theorem,  the  generator  of  the  symmetry  is 


Ek(n)  =^fL^5Uk(n)+T^-5Uk+(n) 

auk(n)  8ljk  (n) 


=  i^-{tr(U^  (n)  ^-Uk(n)  )-h.c.} 

gf  Z  (2-29a) 


and  similarly 


(n)  =  -i  -jL  { tr(Uk  (n)  ^  Uk  (n) )  -  h.c. 

gt  (2-29b) 


By  the  use  of  the  following  identities 


U+  U  =  1  U+U  +  U+U  =  0  (2-30a) 
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and 


N2-l 


I  ^ay<VN5iiVVu 

1=1 


(2-30b) 


it  is  possible  to  compute  the  quadradic  Casmir  operator 


N  -1    a  o  2  + 

£  E  (n)  Eka(n)  =  ^-Ntr(Uk  (n)Uk(n)) 

a=l  gt4 

N2-1  _a  _a 


=    I   E (n)E  (n) 

a=l  K  (2-31) 


In  terms  of  the  Casimir  operator,  the  Hamiltonian  becomes 


2 


H  =  ^-IZEa  E*  +4-I(N-RetrUp) 
aN  !  a  1    1  g2aP 


's  (2-32) 

where  a  constant  has  been  added  to  the  Hamiltonian  such  that  when 

U,  =  1  for  all  k,  H  =  0.  If  the  geometric  mean  of  gs  and  gt  is  defined  as  the  bare 

coupling  constant  gpj  in  Hamiltonian  formulation 


gH    gsgt  (2.33) 
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eq.  (2-32)  becomes  the  Kogut-Susskind  Hamiltonian(5)  of  SU(N)  lattice 
gauge  theory 


J"Iote  that  the  prefactor  —  approaches  unity  in  the  continuous  space  limit, 


so  this  factor  can  be  ignored  for  spectrum  calculations  in  the  continuum. 

2.2.2.  Discussions  and  Comments 

When  the  Kogut-Susskind  Hamiltonian  was  derived  from  the  Wilson 
action  in  the  previous  section,  the  temporal  gauge  first  was  specified  to  start 
with.  But  as  was  mentioned  earlier,  the  Hamiltonian  still  possesses  the 
invariance  of  time-independent  local  gauge  rotation.  This  symmetry  is  related  to 
the  extra  degree  of  freedom  in  gauge  fields.  In  order  to  complete  the  definition 
of  Hamiltonian,  it  is  necessary  to  eliminate  explicitly  the  unphysical  degree  of 
freedom.  However,  in  the  continuous  nonabelian  gauge  theory  there  appear  to 
be  some  special  problems.  These  problems  will  be  discussed  first  in  the  context 
of  continuous  nonabelian  gauge  theory,  then  in  the  context  of  lattice  gauge 
theory. 

In  the  continuous  nonabelian  gauge  theory,  the  classical  Lagrangian  density 

is 


^-I(N-RetrUp)} 


(2-34) 


(2-35a) 
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where 


Ea  =  VA°a-4^-g„fabCAbA0C 

dt     60  (2-35b) 


and 


_a  —  —a 
B  =VxA 


be 


2  too 


-b 

A  xA 


(2-35c) 


The  equations  of  motion  are 


—   -» a         ahc  —  b  _  c 

V.E  +gof*A  •  E  =0  (2.36a) 


and 


^   -  3Fa        abc   bo  -c  _b  -c 
VxB-^-  +  gcf     (A     E+AxB)  =  0 


(2-36b) 


To  cast  the  theory  in  Hamiltonian  form,  notice  that  eq.  (2-36a)  is  the  nonabelian 
version  of  Gauss'  law,  and  A0  is  not  an  independent  variable.  It  can  be 
eliminated  through  eq.  (2-36a)  by  the  use  of  eq.  (2-35b).  Therefore  it  is 
natural  to  choose  the  temporal  gauge  A°  =  0  to  case  the  theory  into  the 
Hamiltonian  form.  The  Hamiltonian  density  is 
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,   -a  _a  _a  _a 


(2-37a) 


with 


(2-37b) 


_a 


and  B  is  the  same  as  eq.  (2-35c).  There  are  two  equations  of  motion.  One 
gives  eq.  (2-37b),  and  the  other  gives  eq.  (2-36b).  However,  Gauss'  law,  eq. 
(2-36a),  is  missing.  Therefore  it  is  necessary  to  impose  Gauss'  law  as  a 
constraint  on  the  Hamiltonian.  In  general,  there  are  two  ways  to  solve  the 
constraint,  and  each  one  has  its  own  difficulty: 

a)  The  constraint  is  transferred  to  the  gauge  condition  to  impose  on  the 
field  variables  before  quantization.  Combining  eqs.  (2-36a)  and 
(2-37b)  together,  Gauss'  law  becomes 


=  0 


(2-38a) 


Since  the  structure  constant  f^30  is  antisymmetric  in  b  and  c 


„abc ?i~a^  i  a  abc-*b  -.b_,c 
gof      A  .M_  =  lgo|(f     A   .A   A  ) 


=  0 


(2-38b) 
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and  Gauss'  law  becomes 


a  —   —  a 
■2-  V»  A  =0 


at 


Therefore,  imposing  the  constraint  is  equivalent  to  working  in  the 
Coulomb  gauge 


—   —  a 

V  •  A   =0.  (2-39) 


However,  Gauss'  law  in  nonabelian  gauge  theory  is  a  nonlinear  equation, 

-a  _a 

solution  to  which  does  not  naturally  separate  E  and  A    into  transverse  and 
longitudinal  parts  to  incorporate  the  Coulomb  gauge.  It  was  also  pointed  out  by 
Gribov(27)  that  there  are  ambiguities  involved  in  extracting  the  dynamical 
variables  in  the  Coulomb  gauge. 

b)  The  constraint  is  imposed  on  physical  states  of  Hilbert  space:  If  a  state 
in  Hilbert  space  is  represented  by  \|/  [A]  which  is  in  coordinate  representation 
then  the  Schroedinger  operator  E  (x)is  represented  by 


E(x)  =  i^A_ 
5A(x) 

Under  the  time-independent  gauge  transformation 


A  V)  ->  A  a(x)  =  A 3  (x)  -  V  I  a  (x)  +  g0 f  abcl  b(x)  A°(x) 


(2-40) 


(2-41) 
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the  wavefunction  changes  to 


Vt  A      -  V  Ea(x)  +  gQf  abc  Lb(x)  A°(x)  ] 


=  exp  {  jd3x  (  V  La(x)  -  g()f abC  Zb(x)  A°(x) )  •  — |—  }  ¥[  A  ] 

5A  (x) 


=  exp  { -iJdS  •Ea(x)Za(x)  +  i  jd3xi:a(x)(5aC  V  +  gQf  abc  a\x)  ) 
•  I  (x)}y[A]  (2_42) 

If  the  surface  term  in  eq.  (2-42)  can  be  thrown  away,  the  constraint  can  be 
satisfied  by  requiring  gauge  invariance  for  physical  states  in  Hilbert  space,  i.e., 

— •  >  — •  ac  >  oKp  »  b  — .  C   , 

ifVlA]=V[A],then{(8    V  +  g0fdDCA   )»E   }  \|/[  A  ]=  0  However, 

because  of  the  non-trivial  topological  structure  of  nonabelian  gauge  fields,  the 
surface  term  cannot  be  thrown  away.  This  leads  to  the  famous  ®  -vacuum  in 
QCD.(28) 

Next  to  be  considered  is  the  situation  in  lattice  gauge  theory.  The  Kogut- 
Susskind  Hamiltonian  is  symmetric  under  time-independent  local  gauge 
rotation.  An  operator  which  performs  such  transformations  can  be  constructed 


Q[Z(n)]  =  exp{iZ££a(n)£(Ea(n)  +  R  (n)  } 

n  a  k     k         k  (2-43) 
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Note  that  by  using  eq.  (2-27)  it  can  be  seen  that,  infinitesimally,  the  operator  £2 

does  generate  the  transformation  of  eq.  (2-26).  The  Hamiltonian,  eq.  (2-34), 

commutes  with  CI.  This  follows  since  the  Casimir  operator, 

a  '  cornmutes  witn     generators  of  the  group.  Secondly,  the  effect 

of  Q  on  each  link      (n)  is  to  perform  the  transformation 

U  (n)  ->  U'(n)  =  exp[i  I  E%)  £  ]  Uu(n)  exp[  -  i  2  Ea  ( n  -  (L )  tr  \ 
*■  y-  a         2      ^  a  2  (2-44) 

so  the  magnetic  term  of  H,  £  (N-RetrUp) ,  is  constructed  to  be  invariant 

P 

under  this  transformation.  Since 


[H,Q]=0  (2.45) 

the  eigenf unction  \|4      ]  is  degenerate  with  the  eigenf unction  ]. 
If  it  is  required  that  the  physical  states  should  be  singlet  under  the 
transformation  of  Q,  i.e., 

nV[uk]  =  v[uk]  (246> 

it  implies  that 


X(E*(n)  +  EJn))y[U  1=0 

k  k        .  (2-47) 
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a  a- 

Using  the  definitions  of  E^(n)  and  (n)  in  eqs.  (2-29a)  and  2-29b),  it  can 
be  shown  that  in  the  limit  a->  0  eq.  (2-47)  becomes 


{(5aCV  +  g0fabcAb)«EC}  y[A]=0  (24g) 


Therefore,  the  constraint  of  Gauss'  law  is  solved  by  requiring  gauge  invariance 
-of  the  wavefunction.  Here  there  is  no  difficulty  caused  by  non-trivial 
topological  structure  of  gauge  fields  because  the  topological  structure  of  gauge 
fields  is  trivial  on  the  lattice,  as  discussed  in  Section  2.1.2. 

Finally,  two  remarks  can  be  made  about  Hamiltonian  formulation: 

(a)  Since  time  is  kept  continuous,  it  is  easier  to  do  numerical  simulations  in 
the  Hamiltonian  formulation  because  the  number  of  degrees  of  freedom  is 
reduced.  But  variables  associated  with  each  link  are  doubled. 

(b)  The  bare  coupling  constant  in  Hamiltonian  formulation  is  defined  as 

gH  =  gS  gt  * 


It  is  not  clear  how  gH  relates  explicitly  to  the  bare  coupling  constant  8L  in 
Lagrangian  formulation.  Since  the  correlation  effect  depends  on  the 
dimensionality  of  the  lattice  and  higher  dimensions  imply  more 
correlations,  it  is  therefore  expected  that  for  the  same  physical  phenomena, 
gH  <  gL'  As  wiU  be  shown  later,  this  is  the  case  for  rough  transition  in 
the  string  tension. 
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2.3.  Continuum  Physics  and  Scaling  Behavior 


As  was  mentioned  often  in  the  previous  sections,  the  continuum  physics 
will  be  recovered  in  the  limit  of  the  lattice  spacing  approaching  zero.* 
However,  this  does  not  mean  that  the  continuum  physics  is  reproduced  by 
shrinking  the  lattice  spacing  to  zero.  This  is  an  essential  point  of  lattice  gauge 
calculation.  When  space-time  is  discretized,  a  cut-off  is  introduced  naturally 
into  the  theory  and  the  cut-off  represents  a  hidden  scale  of  the  theory.  Upon 
renormalization,  this  hidden  scale  will  be  translated  to  an  arbitrary  finite  scale  in 
the  form  of  a  floating  renormalization  point.  If  lattice  spacing  a  keeps  being 
reduced  and  the  coupling  constant  go  keeps  being  renormalized,  different  but 
finite  scales  will  appear.  Since  nonabelian  gauge  theory  is  an  asymptotic  free 
theory,  when  a  is  reduced  so  that  it  is  sufficiently  small  then  g0  is  weak  and  the 
relation  between  gG  and  a  is  calculable. 

The  (3-  function  in  Wilson's  action  can  be  defined  as 


p  =  a 


dg0(a) 


da 


(2-49) 


In  weak  region,  the  perturbative  calculation  gives 


P(g0)  =  |30g03  +  P1g05  +  O(g7) 


(2-50a) 


♦This  section  summarizes  and  somewhat  expands  a  discussion  found  in 
ref.  (30),  it  includes  a  number  of  results  which  will  be  used  in  the  later 
chapters. 
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where  (3Q  for  SU(N)  gauge  theory  without  fermions  is  given  U)  by 


Po  =  (_l^)(IiN) 

16*        5  (2-50b) 


to  one-loop  correction,  and  Pj  is  from  a  two-loop  correction  (29)  with 


P^-^Wn2) 
1      16kz       *        .  (2-50c) 


In  general,  (J (go)  depends  on  the  renormalization  scheme.  However,  the  first 
two  coefficients  in  the  asymptotic  series,  PQ  and  P  j,  are  universal.  Combining 

eqs.  (2-49)  and  (2-50)  together,  we  have 


=  d(log  a) 


P0So  +  Pi  4  +  °(go) 

Integrating  both  sides  gives 

Pi 


(2-5 la) 


ATa  =  (^T)2P°exp{-^-} 
Po4  2P^ 


(2-5  lb) 


where  is  a  constant  of  integration.  Some  remarks  may  be  made  about  this 
important  result,  eq.  (2-5 lb): 
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(a)  The  result  is  up  to  the  two-loop  correction,  and  is  independent  of  the 
renormalization  scheme,  because  PQ  and  p  ^  are  independent  of  the  scheme. 

(b)  In  order  to  define      properly,  a  two-loop  correction  is  needed.  The 

two-loop  correction  gives  the  factor  before  the  exponential  decrease  of  the 

coupling  on  the  right  hand  side  of  eq.  (2-5  lb)  and  in  turn  gives  a  proper  and 
finite  value  of  A^. 

(c)  The  value  of     is  scheme-dependent  and  the  relation  between  two 

's  can  be  obtained  through  a  one-loop  calculation  in  perturbation  theory.(30) 
In  order  to  show  this,  consider  two  difference  schemes  each  defining  a  bare 
coupling  as  a  function  of  cut-off:  gQ(a)  and  gc  (a).  They  must  agree  with  each 
other  in  the  classical  limit 


g0  =gG  +  cgo 


(2-52) 


so  then 


Pi 


1 


Pi 


2Po4 


1 


}  exp{  } 


o 


(2-53) 
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Therefore 


a't 

AL  Po 


(2-54) 


Hasenfratz  and  Hasenfratz(31)  performed  the  necessary  one-loop  calculation 
and  found 


a  .  83.5  For  SU(3) 

_R  < 

A  1  57.5  For  SU(2) 

^  (2-55) 


and 


0.91    For  SU(3) 


AT       '0.84  ForSU(2) 
L  (2-56) 


where  A  is  calculated  in  the  Feynman  gauge  for  continuous  gauge  theory,  and 

A„  is  in  the  Kogut-Susskind  Hamiltonian. 
rl 

Now  consider  that  a  physical  quantity,  glueball  mass  M  for  example,  is 

s 

calculated  on  a  lattice.  If  at  a  certain  range  of  coupling  constant  it  is  observed 
that  the  dimensionless  quantity  M  a  scales  with  respect  to  the  coupling  g0  as 
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A. 

/  1  N  .  1  , 

M    a  =  constant  •  ( — — — )       exp  {  } 


then  the  correct  continuum  physics  has  been  reproduced  and  an  intrinsic  scale  is 
set  in.  In  the  example  here,  the  constant  gives  the  glueball  mass  in  the  unit  of 
intrinsic  scale  A    Such  a  range  of  go  is  called  the  scaling  window. 


CHAPTER  3 
ANALYTICAL  CALCULATION  METHODS 


This  chapter  is  devoted  to  discussions  of  analytical  methods  used  in 
lattice  gauge  theories  and  to  establishment  of  a  systematic  way  of  calculating 
the  physical  observables.  As  mentioned  in  Chapter  1,  the  analytical 
"calculation  is  usually  done  by  using  strong-coupling  expansion  both  in 
Lagrangian  and  Hamiltonian  formalisms.  In  Section  A,  the  calculation  of 
string  tension  is  used  as  an  example  to  discuss  the  strong-coupling  expansion. 
In  the  next  two  chapters,  the  variational  method  is  used  to  perform  practical 
calculations  in  Hamiltonian  formalism.  To  set  the  stage  for  later  calculations, 
a  systematic  way  of  computing  the  physical  observables  analytically  in 
variational  calculations  is  presented  in  Chapter  2.  The  method  is  a 
combination  of  Robson's  plaquette-space  integration^2)  and  Batrouni's  mean 
plaquette  metho&O 1)  Although  the  discussions  of  Section  B  are 
concentrated  to  the  Hamiltonian  formalism,  the  method  also  is  applicable  to 
calculations  in  Lagrangian  formalism 

3.1.  Strong-Coupling  Expansion 
The  strong-coupling  expansion  is  a  common  and  conceptually  straight- 
forward way  of  performing  analytical  calculations  in  both  Lagrangian  and 
Hamiltonian  formalisms,  and  the  results  for  various  physical  observables  have 
been  obtained  by  many  authors.(22>27>4°)  As  mentioned  in  Chapter  1,  on  a 
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large  lattice  spacing  the  theory  is  in  a  strong-coupling  regime,  and  the  strong- 
coupling  expension  technique  can  be  employed.  In  this  section,  the  expansion  of 
the  string  tension  in  the  order  of  the  coupling  strength  is  used  as  an  example  to 
discuss  the  method.  The  calculations  presented  here  are  only  carried  out  to  the 
non-trivial  leading  order  in  both  the  Lagrangian  and  the  Hamiltonian  formalisms, 
and  the  results  are  somewhat  well  known.  The  purposes  of  presenting  these 
simple  calculations  are  twofold:  One  is  to  show  how  the  linear  confinement  arises 
naturally  in  the  lattice  gauge  theories.  Particularly,  the  physical  picture  of  linear 
confinement  is  clearly  established  in  the  non-trivial  leading  order  calculation  of 
strong-coupling  expansion  of  the  static  quark  potential  in  the  Hamiltonian 
formalism.  The  physical  picture  thus  obtained  is  helpful  for  the  calculations 
presented  in  Chapter  4.  Secondly,  through  the  calculation  of  the  string  tension, 
the  limitation  of  the  applications  of  strong  coupling  expansion  can  be  discussed 
clearly. 

In  Lagrangian  formalism,  the  strong  coupling  expansion  is  a  perturbative 
expansion  in  1/g^  .  To  obtain  the  strong-coupling  expansion,  consider  the 

expectation  value  of  a  rectangular  Wilson  loop  of  size  I  •  J  (Figure  3)  in  SU(N) 
lattice  gauge  theory 


<  W(C)  >  = 


[ [ {?dUk]  (  rT  Vfec  Vn)  >CXP^f  (  *UP  +h-C-  j ] 
J[ndUk]exp[^2(trUp  +  h.c. )] 


(3-1) 
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Here  the  normalization  constant  in  the  action  has  been  cancelled  between  the 

numerator  and  the  denominator,  and  the  factor     is  inserted  into  the  Wilson  loop 

for  convenience.  Not  that  the  Haar  measure  ndU,  is  normalized  such  that 

k  k 


ndU,  =1 

k     k  (3-2) 


-Since 


[IldU,  ]UL,(n)  =  0 
k     k     ^  (3-3) 


it  is  observed  that  eq.  (3-1)  vanishes  at  the  zero111  order  in  p\  To  go  to  the  first 
order  in  p\  every  link  from  the  action  must  have  a  partner,  either  from  the  action 

itself  or  the  Wilson  loop.  To  fill  out  the  Wilson  loop  with  a  particular  orientation 

(Figure  4),  we  should  have  I  •  J  plaquettes  from  the  action,  which  then  gives  the 
/  p  vI«J 

factor  ^  — - )         .  For  any  two  overlapped  links  in  opposite  directions 
JdUU.-U.'^-i-S..  5., 

J        ij    kl     N  il  jk  (3_4) 
where  i,  j,  k,  and  1  are  indices  of  matrix  element.  All  together,  there  are  [21  •  J  + 


(I  +  J)]  links  within  the  Wilson  loop,  including  those  in  the  boundary.  The  link- 

r     „  ,x   •        ,       /  1  x  21  •  J  +  (I  +  J) 
integraaon  of  eq.(3-4)  gives  a  factor       )  .  After  the  link 

integration,  the  group  element  associated  with  each  link  becomes  an  identity,  and 
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each  trace  gives  a  factor  N.  In  total,  there  are  [  I  •  J  +  (I  +  J)  +  1]  traces  within 
the  Wilson  loop  including  the  boundary.  Therefore,  there  is  another  factor, 
n!*J  +  (I  +  J)+1    Combining  all  of  the  factors  together,  the  result  is 


<wrn>  -  ±-t  p  \w<  1  .2I-J+a+J)  l-J+(l+J)+l 


2N2 


for  SU(N),  N>2. 


(3-5a) 


t 

J 

I 


Figure  3.  A  Wilson  loop  of  size  I  •  J. 

Since  SU(2)  representation  is  real,  it  is  not  necessary  to  be  concerned  about  the 
orientation.  The  factor  v  2N  should  be  replaced  by  v  N  '       ,  and 

the  result  becomes 


<wroVk)<7)W 


for  SU(2). 


(3-5b) 
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I 


Figure  4.  A  Wilson  loop  filled  out  by  plaquettes  with  proper  orientations. 


At  this  non-trivial  order  of  strong-coupling  expansion,  an  area  dependence  of  the 
Wilson  loop  has  been  seen  already 


<W(c)>  >  ~  exp  [  -a  area  (c)  ] 


(3-6) 
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where  a  is  string  tension  and 


log(-^) 


for  SU(N),  N>2 


2N' 


a  a 


(3-7) 


for  SU(2) 


The  procedure  can  be  continued  to  obtain  high  order  terms  in  the  expansion  and 
jhen  try  to  extract  the  continuum  physics  from  the  expansion.  Before  the 
discussion  of  continuation,  it  is  important  to  observe  what  happens  in 
Hamiltonian  formalism. 

The  quantity  in  Hamiltonian  formalism,  related  to  the  Wilson  loop  in 
Lagrangian  formalism,  is  the  static-quark  potential.  The  area  dependence  of 
Wilson  loop  implies  that  the  static-quark  potential  is  linear  to  the  distance  between 
a  pair  of  quark  and  antiquark.  Therefore,  the  static-quark  potential  between  a 
heavy  quark  at  site  n  and  a  heavy  antiquark  at  site  n  +  R  is  considered.  To 
simplify  the  calculation,  the  gauge  group  is  specified  as  SU(2).  The  Kogut- 
Susskind  Hamiltonian  of  SU(2)  lattice  gauge  theory  is 


H  = 


(3-8a) 


As  mentioned  in  Chapter  2,  the  Hamiltonian  formalism  also  permits  a  strong- 
coupling  expansion  which  is  now  an  application  of  quantum  mechanical 
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perturbation  theory.  In  the  strong-coupling  limit,  the  Hamiltonian  of  eq.(3-8)  is 
dominated  by  the  color-electric  term 


g2 

Ho  =  ^?iEkEk 

*  ka  K  K.  (3-8b) 


For  each  link,  it  is  a  quadratic  Casimir  operator.  The  corresponding  vacuum  of 
eq.  (3-8b)  thus  has  each  link  in  a  color  singlet  state 


Ef  I  0  >  =  0 

k  .  (3-9) 


Now  consider  a  heavy  q  q  state,  a  gauge  invariant  operator  is  needed  to  describe 
the  state.  The  simplest  choice  of  the  operator  is 


¥(n+R)[    n    U.  ]¥(n)IO> 

kepath   k  (3_1Q) 


where  the  path  extends  from  n  to  n  +  R  with  the  minimum  number  of  links 
required,  *F  (  n )  is  the  operator  which  creats  a  quark  q(n)  at  site  n,  and 

(  n+R )  creats  an  antiquark  q(n  +  R)  at  site  n  +  R.  Note  that  the  color  indices 
are  contracted  and  supressed.  For  a  link  there  is 
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E E?  E.a U.  IO>  =  (Z^-»^)U,  I  0>  =  4-U.  I0> 

akkk  va2      2      k  4k  (3-11) 


where  x  a,  a  =  1,2  and  3  ,  are  the  Pauli  matrices.  There  is,  correspondingly,  the 
potential  between  a  quark  and  an  antiquark  separated  by  distance  R, 


V(R)=Si«l«^=oR 

2a    4    a  (3.12) 


where  y  is  the  minimum  number  of  links  required  in  the  path  of  eq.  (3-10).  The 
string  tension  then  is 


g?,»i8  H 

H  .  (3-13) 


The  result  obtained  here  is  the  leading  term  in  the  strong-coupling  expansion.  In 

principle,  quantum  mechanical  perturbation  theory  can  be  applied  to  the 

2 

calculation  to  obtain  higher  order  terms  in  the  power  series  of  1/  As 
stressed  in  Chapter  2,  to  extract  the  continuum  physics,  the  physical  observables 
calculated  should  exhibit  scaling  behavior  in  a  certain  range  of  coupling  constant. 
Since  the  scaling  usually  happens  in  the  intermediate  coupling  region,  it  may  be 
expected  that  if  the  higher  order  terms  in  the  expansion  are  included  the  scaling 
behavior  will  appear.  In  the  case  of  string  tension,  for  large  enough  M  it  is 
possible  to  extract  the  continuum  physics  from  the  expansion 
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(3-14) 


where  P  =  ~L ,  and  fQ  (     )  is  equal  to  log  ( )  in  Lagrangian  formalism  and 
go 

^-  g^  in  Hamiltonian  formalism.  However,  in  the  following  two  situations  the 
8  H 

jstrong-coupling  expansion  is  not  sufficient  to  determine  the  continuum  physics: 

(a)  If  any  critical  point  exists  between  the  strong-coupling  and  the  scaling 
regions,  the  expansion  coefficients  will  fluctuate  largely,  and  it  is  necessary  to 
continue  beyond  the  critical  point  in  the  expansion  in  order  to  obtain  the  seating. 

(b)  If  the  scaling  region  happens  at  the  weak-coupling,  or  if  the  symmetry 
(used  to  assign  the  quantum  state  of  the  calculated  physical  observable)  can  only 
be  restored  in  the  weak-coupling,  the  expansion  will  not  be  able  to  give  the 
continuum  physics. 

To  save  the  strong-coupling  expansion  in  these  situations,  it  is  possible  to  employ 
the  Pade  approximation  (9)  when  continuing  the  result  of  strong-coupling 
expansion  beyond  the  critical  point  or  going  into  the  weak-coupling  region. 

The  Pade  approximation  has  been  studied  in  statistical  mechanics  for  a  long 
time.  The  main  idea  of  this  approximation  is  to  replace  the  power  series  in  the 
expansion  by  the  ratio  of  two  power  series 


(3- 15a) 
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here  N  +  D  ^  M.  The  coeffiecients  ni  and  di  are  determined  by  solving  the 
equation 

D      i      M        n      N  i 
(Id  p  )(  I  bnP    )  =  I  ni 

i=o         n-o  j=o  J  (3_15b) 

"order  by  order.  The  critical  point  corresponds  to  pc  with 


D  „i 
£  &p  =0 

i=0  1  c      .  (3-15c) 


However,  if  the  system  exhibits  an  unphysical  singularity  with  convergence 
radius  larger  than  that  of  the  physical  one,  the  straightforward  application  of  the 
Pade  approximation  can  not  give  the  correct  information.  It  is  necessary  to  make 
a  transformation  of  the  variable  P  such  that  the  convergence  radius  is  determined 
by  the  physical  singularity.  Unfortunately,  this  is  also  the  case  for  the  string 
tension  in  Lagrangian  formalism  and  it  is  not  known  how  to  transform  away  the 
unphysical  singularity  which  is  the  rough  transition  point.(33) 

From  the  leading  order  calculation  of  the  static  quark  potential,  it  can  be  seen 
that  in  the  strong-coupling  region  the  formation  of  color-electric  flux  tubes  gives 
rise  to  a  linear  potential  between  the  quarks  and  supports  the  confinement 
phenomena.  The  quantum  effect  will  set  fluctuatation  in  the  flux  tube  and  the 
transversed  width  of  the  flux  tube  is  inversely  proportional  to  the  glueball  mass. 
The  rough  transition  of  the  string  tension  corresponds  to  the  point  where  the  flux 
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tube  starts  to  dislocate  and  continuous  translation  symmetry  is  restored 
effectively.  The  rough  transition  here  is  not  a  physical  singularity  because  the 
shape  of  the  flux  tube  is  still  maintained  even  past  the  transition  point. 

In  the  Lagrangian  formalism,  Monte  Carlo  results  show  that  the  rough 
transition  point  exists  between  the  strong-coupling  and  the  seating  regions  for 
SU(2)  and  SU(3)  (4)  and,  as  mentioned  in  the  beginning  of  this  section,  it 
raises  serious  questions  about  the  continuation  of  strong-coupling  expansion 
beyond  the  rough  transition  point.(29)  Our  calculation  results  of  Chapter  4 
show  that  a  similar  situation  also  appears  in  the  Hamiltonian  formalism 
Therefore,  it  also  is  questionable  to  use  the  strong-coupling  expansion  to 
calculate  the  string  tension  in  the  Hamiltonian  formalism.  Considerable  effort 
have  been  concentrated  recently  on  variational  calculations  which  are  able  to 
circumvent  these  difficulties,  since  they  find  the  estimates  for  the  physical 
observables  on  both  sides  of  these  "phase  changes." 

3.2.  Plaquette-Space  Integration  and  Mean  -  Plaquette  Approximation 
This  section  describes  the  analytical  method  used  to  evaluate  the  integrals 
necessary  for  variational  calculations.  The  idea  behind  variation  methods  is 
simple:  A  trial  wavefunction  with  variational  parameters  is  proposed  for  a 
physical  state  and  the  parameters  then  are  fixed  by  mmimizing  the 
corresponding  energy.  The  calculations  that  must  be  carried  out  are 
computations  of  expectation  values  of  physical  observables  with  respect  to  the 
trial  wavefunction.  Either  the  Monte  Carlo  method  or  the  analytical  method  may 
be  used  to  perform  the  calculations.  A  systematic  way  of  performing  the 
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analytical  calculations  will  be  developed  in  the  next  section  and  use  of  the  Monte 
Carlo  method  in  solving  the  problem  will  be  the  subject  of  Chapter  6. 
3.2.1.  Method 

In  the  Lagrangian  formalism  the  Wilson  action  is 


S(U_)  =  X^(trUD  +  h.c.) 

V     P  g        P  .  (3-16) 


The  ground  state  expectation  value  of  a  physical  observable  O  ( Up  )  is  given 


by 


<0(UD)>  =  lJ  [DU]0(U_)exp[-S(U_)] 

v        L.  v  v  (3- 16a) 


where  Z  is  the  vacuum  to  vacuum  transition  amplitude 

Z  =  J[DU]exp[-S(Un)] 

V  P  (3-16b) 

and  DU  is  the  Haar  measure  defined  in  link  space 


DU  =  n  dU, 
k  k 


which  satisfies  the  left-  and  right-invariance 


(3- 16c) 


DU  =  D(VU)=D(UV) 


(3-16d) 
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where  V  is  an  arbitrary  element  of  the  SU(N)  gauge  group.  Because  the 

physical  observable,    ^p ' ,  is  a  function  of  the  traces  of  the  various  kinds  of 

plaquettes,  exists,  therefore,  is  a  functional  integral  in  which  the  integrand  is  a 

function  of  the  plaquette  and  the  measure,  described  in  eq.  (3- 16c),  is  defined  in 

link  space.  In  1979,  Halpern(34)  found  a  unique  inversion  for  the  potential 
(x)  in  terms  of  the  field  strength  F^  (x)  in  continuum  gauge  theories  by 

using  completely  fixed  axial-like  gauge.  The  variables  of  the  Haar  measure 

therefore  can  be  changed  from      (x)  to       (x).  After  fixing  the  gauge 

completely,  A^  (x)  has  two  degrees  of  freedom,       (x)  has  six  degrees  of 

freedom.  The  constraint  is  consequently  required  to  be  imposed  on  the  field 
configurations  of  F^  (x).  Halpern(34)  showed  that  the  constraint  is  the 

Bianchi  identities 


(3-17) 


where  p,,   =  J-G         F^^  is  the  dual  field  strength.  The  Bianchi  identities 

have  four  components  at  each  space-time  point.  The  constraint  therefore 
reduces  the  degrees  of  freedom  of  F^  (x)  to  2.  However,  not  all  of  the 
components  of  the  Bianchi  identities  are  independent  of  each  other.  Consider  a 
completely  fixed  axial-like  gauge  such  as 


A3  (x,y,z,t)  =  A^x.y^)  =  A2  (t,x0,y,z0)  =  AQ  (t,x0,y0,z0)  =  0 
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where  ^o'^o^o^o^  is  a  particular  space-time  point  chosen  to  start  the  gauge 
path.  The  third  component  of  Bianchi  identity  is  only  required  at  z  =  zq 

*3  (x'y' V>  =  °-  (3.19) 
Halpern(31)  called  this  the  3  •  1  Bianchi  identities.  In  general, 


J[DA^]8(CFG)f(F^) 


=  J[DF^]8[I(F^)]f(F^)  (32Q) 

where  5(  CFG  )  is  the  delta-function  which  imposes  the  gauge  condition, 
eq.  (3-18),  and  5[  I(  F^)  ]  is  the  delta-function  for  imposing  the  Bianchi 

identities  to  the  various  field  configurations  in  field- strength  space. 

Following  Halpern's  work,  BatrouniO  1)  gave  the  lattice  version  of  the 

Bianchi  identities  as 


6[I(F^)]^o(Pc-l) 


By  making  the  gauge  transformation  Up  ->  Up,  the  lattice  Bianchi  identity  can 


be  abelianized  as 


5(Pc-l)->5(Pc-l)  =  8(nU  (i)-l) 

i=l  (3-22) 
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where  U    ( 1 ) ,  •  •  •  • ,  Up  (6)  are  the  six  plaquettes  which  form  a  cube  C  and  the 
orientation  is  assumed  to  be  taken  into  account  properly.  However,  the  gauge 
degrees  of  freedom  are  completely  fixed  during  the  course  of  the  unique 
inversion  from      to  Up  .  The  question  of  abelianizing  the  Bianchi  identity  is  a 
question  of  the  unique  correspondence  between  a  gauge-fixed  plaquette  and  a 
gauge- invariant  plaquette.  By  constructing  explicitly,  Batrouni(l  1)  showed  that 
the  abelianization  can  be  done  over  large  subspaces  of  the  lattice  but  will  not  be 
complete. 

The  above  arguments  can  be  applied  to  the  Hamiltonian  formalism  of  lattice 

gauge  theories  equally  well.  In  the  Hamiltonian  formalism,  the 

Kogut-Susskind  Hamiltonian  is  obtained  from  the  temporal  gauge,  and  the 

physical  states  are  subjected  to  the  nonabelian  Gauss'  law,  which  implies  that 

the  physical  states  are  local  gauge  invariant.  In  the  coordinate  representation, 

the  wavefunction  is  a  function  of  the  traces  of  the  various  plaquettes.  For  a 

local  gauge  invariant  observable  0(Up),  the  expectation  value  with  respect  to 
the  wavefunction  V(U  )  is 


<0(Up)>  = 


J  [DU]  ¥+  (Up)  0(Up)  X|/(Up) 
J[DU]V+(Up)v(Up) 


(3-23) 


which  is  quite  similar  to  eq.  (3- 16a).  To  perform  the  integration  of  eq.  (3-23) 
in  plaquette  space,  it  is  possible  to  choose  the  completely  fixed  temporal-like 
gauge 
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AQ(t,x,y,z)  =  A3(tQ,x,y,z)  =  A2(t0,x0,y,z0)  =  A1(t0,x,y,z())  =  0 


24) 


The  temporal-like  gauge  is  compatible  with  the  temporal  gauge  Aq  (t,x,y,z)  =  0, 
which  is  the  gauge  used  to  obtain  the  Kogut-Susskind  Hamiltonian.  Similar  to 
the  Lagrangian  formahsm,  the  lattice  Bianchi  identities  must  be  imposed,  and 
the  zerom  component  of  the  Bianchi  identities  is  only  required  at  t  =  to 

ll  Mo)  =  0-  (3_25) 

The  difference  is  in  the  number  of  components  of  the  Bianchi  identities  that 
must  be  imposed  for  a  particular  plaquette.  In  the  Hamiltonian  formalism,  the 
Haar  measure  is  defined  in  three  spatial  dimensions.  When  the  measure  defined 
in  the  link  space  is  converted  to  that  in  the  plaquette  space,  the  Jacobian  is  the 
three  spatial  components  of  the  Bianchi  identities.  This  is  different  from  the 
Lagrangian  formalism,  where  the  Jacobian  is  equal  to  the  four  components  of 
the  Bianchi  identities. 

Working  in  plaquette  space,  as  pointed  out  by  Batrouni,(l  1)  the  Bianchi 
identities  can  be  seen  as  the  gradual  restoration  of  strong  coupling  expansions. 
For  the  nonabelian  lattice  gauge  system,  it  has  been  demonstrated  by  various 
numerical  calculations  that  the  confined  phase  coexists  with  the  non-confined 
phase,  and  that  the  only  fixed  point  is  at  g  =0.(4)  Moreover,  it  has  been  shown 
and  argued  that  the  vacuum  magnetic  fluctuations  are  quite  localized.^) 
Because  of  the  localization  of  the  magnetic  fluctuations  in  the  vacuum,  it  is 
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justifiable  to  work  in  a  local  region  with  the  abelianized  lattice  Bianchi  identities 
while  neglecting  the  nonabelianized  identities  which  usually  occur  in  several 
lattice  spacing  units  away  from  the  local  region.  Using  the  property  of 
localization,  Robson(12)  Was  the  first  to  combine  the  Bianchi  identities  with 
character  expansion  in  group  representation  to  establish  an  analytical  calculation 
method.  The  essential  point  of  the  method  is  the  need  for  using  a  small  number 
of  the  Bianchi  identities  to  obtain  the  accuracy  required.  In  the  following 
section,  Robson's  methodO^)  js  modified  by  combining  it  with  Batrouni's 
mean-plaquette  method(^)  and  the  method  is  examined  in  detail  for  the  case  of 
the  SU(2)  gauge  group.  The  method  then  may  be  applied  in  a  straightforward 
manner  to  the  other  gauge  group. 

The  group  structure  of  SU(2)  is  simple.  Since  the  trace  of  a  plaquette, 
trUp,  is  taken  in  the  fundamental  representation,  Up  is  a  two  by  two  matrix  and 
can  be  parameterized  as 


U-,  =  cos(  -§-)  -  i  n*x  sin( ) 
P  2  2  (3-26) 


where  t*  are  the  Pauli  matrixes  and  na  =(sin8  cos<|),  sinB  sin<|>,  cos9)  is  a  three 
dimensional  unit  vector  with  0  <  p  <  2k,  0  <  8  <  2k,  and  0  <  <j)  <  2k  .  The 
character  of  jth  representation  is  denoted  as  Xj  (Up)  with 
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with  Xj(Up)  satisfying  the  orthonormal  condition 

WW^  (3.27b) 
and  the  product  of  Xj(Up)   can  be  written  as  a  Clebsch-Gordon  series 

m=P"Kl  ,  (3_27c) 

From  eqs.  (3-26)  and  (3-27a), 


trUp  =  x1(Up)  =  2cos(|-) 


To  simplify  the  calculation,  the  product  of  single  plaquette  wavefunction  is 
used  as  the  trial  ground  state  wavefunction 

O  =  exp[  XL  trU  ]  =  n  exp[  torU_  ] 

P     P      P  P  (3_28) 

where  X  is  a  variation  parameter.  The  trial  ground  state  wavefunction  is  exact  in 
the  strong  coupling  limit  and  has  been  shown  to  be  reasonable  in  the  scaling 
region  because  of  the  localization  of  the  magnetic  fluctuations  in  the  vacuum 
Since  G>0  is  invariant  under  the  transformation 


up  ->  vupv+ 


(3-29) 
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where  V  is  an  element  of  SU(2)  gauge  group,  the  abelianized  Bianchi  identity 
can  be  written  as  a  series  of  separable  products 


8(PC-1)=    I   L_6x;(Un) 

T 


'j=0,i  1,-     (2j+l)4i=1"r  Pi 


=  1+1   l——r  [fi  %.  (  U„  )] 

-  -  -  ( 2j+l  )4  pi 
2      2  (3-30) 


In  principle,  the  Bianchi  identities  can  be  expanded  and  always  yield  a  series  of 
various  integrations  effectively  in  two  dimensions.  In  two  dimensions  with  free 
boundary  condition,  the  Jacobian  of  transforming  the  measure  from  link  space 
to  plaquette  space  is  identity  and  the  integration  becomes  trivial.  For  the 
expectation  of  plaquette, 


J  dUn  ( trUp )  exp(  2  X  trUp )  j 


<tp>S<trUp>  =  — f 

J  < 


dU_  exp(  2  X  trU_  )  X  L(4A) 

P  P  r  (3-31a) 


where  ln  (z)  is  the  modified  Bessel  function  defined  as 


In  (z)  =  ^     d0  cosn8«exp(  zcos9  ) 


(3-3  lb) 


57 


and  n  is  an  integer.  By  continuously  differentiating  eq.  (3-3 la)  with  respect  to 
2X,  the  expectation  values  of  various  powers  of  trUp  can  be  obtained.  The 
results  are  listed  in  Table  1.  To  obtain  the  expectation  value  of  a  plaquette  in 
three  dimensions,  the  Bianchi  identities  must  be  inserted  to  change  the  measure 
from  link  space  to  plaquette  space.  In  the  case  of  two  Bianchi  identities  (see 
Figure  5a), 


J[.ftdUn]  [?Z  l— 7  ^ 

P 


i=1    Pi    Jij2(2j1+1)4  (2j2+l)4 


J[.n 


dUn  ]  [II 


i-i    Pi    jiJ2(2j1+l)4  (2j2+l)4 
X-  (U  )X:  (Up  )  {ft  x-  (U   ) }  (Ax.  (U   )  }  ] 

Jt     ^1    J2    n    k=2  Jt    \       n-l  J2  Pn 

X-  (U    )X:  (Up  )  {  ft  X:  (Up  )  }  fllY:  (Up  )  }  ] 
Jl     n    J2    ^1    k=2  Jj     Pfc       n-7   J2  Pn 


[trU  ]  ji  exp(2XtrU  ) 
Pi  m-i  pm 


II  exp(2XtrUn  ) 

P™  (3-32a) 
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2^7 


(a) 


(b) 


Figure  5.  Aplaquette 

a)  with  two  Bianchi  identities;  b)  with  12  Bianchi 
identities. 


If  the  numerator  and  the  denominator  of  the  right  hand  side  of  eq.  (2-21)  are 
both  divided  by  a  common  factor  F 


FsJrJldU D]  fl  exp(  2XtrUt 
1=1    Pj  m=i    r  f 


) 


m 


(3-32b) 


then  both  the  numerator  and  the  denominator  become  a  power  series  of  the 
values  in  Table  1.  For  a  fixed  number  of  Bianchi  identities,  the  convergence  in 
the  expansion  in  j  is  obtained  easily.  In  the  truncation  in  j  the  expectation  values 
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2 

of  [  Xj  ( Up  )  /  ( 2j+l  )3  ]  in  two  dimensions  with  j  =     1     2  ,  •  •  •,  are 
measures  of  the  convergence.  This  quantity  vs.  the  variational  parameter  X  is 
plotted  in  Figure  6. 
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Table  1.  The  expectation  values  of  various  powers  of  trUp  in  two 
dimensions. 


Physical 
Observables 

Expectation  Values  (  t  =  4A.  ) 

trUp^Up) 
2 

t  Ij(t) 

2 

m,24^    12  !0(t) 
(      t2>"  <  V0 

t3up  =  z](up) 

2 

(40  ,  192       g    96.  !0(t) 
■(t  +  t3  )  +  (8+|2)'l1(t)' 

■r^p-xtOJp) 
2 

(16+  432  +  1920  }.(96  +  960  }  !0(t) 
t2       t4         1       t3  V° 

2 

,256    5376    23040  v 
t        t3  t5 

,  ™  ,  1248  ,  11520  ,  !0(t) 

(  ?  t4  v 

2 

(61  ,  4032    76800  322560x 
t2         t4  t6 

,576     18240     161280  x  !0(t) 
t         ■>          t5  jIl(t) 

7  7 

.  .  »  T  T                . ,  /     /"IT  \ 

ir  Up  =  X1(Up) 
2 

,1408  ,  68736  ,  1244160  ,  5160960  * 

"(  t     ''     3            5             7  } 
r         r  t 

KP8  +  10368  i  299520  |  2580480,  r0(t) 
t2        t4         t6  vo 
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Table  2.  The  expectation  value  of  a  plaquette  with  two  Bianchi  identities 
Notations: 

^^VlD'  bs<tr2V2D'  CS<tr3V2D 
d^<tr4Up>2D,  e^<tAjp>2D 

NjSN(I  0)  +  N(0,l)  +  N(I  1) 
2 

=  la4b+JLa9c 
8  256 

Nj  =  N(l,  0)  +  N(0,  1)  +  N(l,  1)  +  N(l  1)  +  N(l,l) 
=  ^(c-a)(b-l)5+^(d.b)a4(b-l)5 

+  65kr(b-1)10(e-2c  +  a) 


Table  2.  Continued 


DjSD(l  0)  +  D(0,l)  +  D(i  1) 


=  la6  +  UJ0b 


8  256 


Dj  =  D(l,  0)  +  D(0,  1)  +  D(l,  I)  +  D(i  1)  +  D(l,  1) 


i(b-1)6+6i8a5(b-1)5(c"a) 


+  65k^(e-2c  +  a)(b-1)10 


If  both  Bianchi  identities  truncate  atj=y,  N  =  N^,D  =  D^. 


If  both  Bianchi  identities  truncate  at  j  =  1,  N  =  N,  +N.,D  =  D1  +  D, . 

J  11  11 
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In  the  figure  it  is  shown  that  the  convergence  in  j  is  very  fast  for  small  X  (strong 
coupling  region)  and  slows  down  for  large  X  (weak  coupling  region).  To 
signify  the  effect  of  the  Bianchi  identities,  for  any  physical  observables  O(Up)  a 
quantity  ,  B,  is  defined  as 

<0<Up)>3D  =  <0(Up»2D(1+B)  (3Ba) 


so 


<Q(Up)>3D-<o(up)>2D 

<0(Un  )>„_ 

P    2D  (3.33b) 


where  the  subscripts  3D  and  2D  denote  that  the  expectation  value  is  evaluated  in 
three  and  two  dimensions  respectively.  Since  the  quantity  B  solely  represents  the 
effects  of  the  Bianchi  identities,  the  analytical  calculation  result  of  the  expectation 
value  of  any  physical  observable,  0(Up),  can  approach  to  the  Monte  Carlo 
calculation  result  when  the  number  of  Bianchi  identities  is  increased  in  the  analytical 
calculation  until  the  convergence  appears  in  the  curve  of  B  vs.  X.  The  analytical 
expressions  of  the  expectation  value  of  a  plaquette  containing  two  Bianchi  identities 
for  the  various  truncations  in  j  is  shown  in  Table  2.  From  the  curves  of  B  vs.  X  in 
Figure  7,  one  can  see  the  result  of  two  Bianchi  identities  converging  at  j  =  1  for  X 
ranged  from  0.05  to  1.1  (or— i- ranged  from  0.1  to  0.5). 
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The  next  step  is  to  increase  the  number  of  Bianchi  identities  up  to  12, 
where  the  B  vs.  X  curve  also  converges  at  j  =  1  for  the  same  range  of  X  as 
shown  in  Figure  8.  A  comparison  between  the  B-curves  of  two  and  12 
Bianchi  identities  is  shown  in  Figure  9,  where  there  is  a  non-negligible 
difference  between  the  two  approximations. 

To  implement  this  method  and  more  rapidly  achieve  the  convergence,  the 
method  can  be  combined  with  Batrouni's  mean-plaquette  method.  By  replacing 
the  surface  plaquettes  of  the  Bianchi  identities  with  the  expectation  value  of 
these  plaquettes, 


V«VUfl^0P>-Vi<rt,P>)  (3-34, 


where  ^j^^V  is  the  Che,Jyshev  polynomial  of  the  second  kind.(36) 
Similar  to  the  previous  procedure,  the  integration  is  then  done  over  the  internal 
plaquettes.  Note  that  in  order  to  obtain  a  consistent  result,  the  number  of  the 
Bianchi  identities  should  not  be  too  small.  In  the  following,  we  will  consider 
two  cases:  One  is  the  expectation  value  of  plaquette  and  the  other  is  the 
expectation  value  of  extended  plaquette. 

(a)  The  expectation  value  of  plaquette:  Twelve  Bianchi  identities  are  used  to 
find  the  expectation  value  of  plaquette  (see  Figure  5b).  After  replacing  the  surface 
plaquettes  by  their  expectation  values  and  carrying  out  the  integrations  for  the  internal 
plaquettes,  it  yields  a  polynomial  equation 
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P(X)  = 


hh  (2j1+ir  (2j2+D 


^MUjjy^i]  Q4[jrj2^,P(X)]  RUrj2,X,P(X)] 


hh  (2j1+l)H  (2j2+l) 


with  P(X,)  =  <  trUp  >,  and  where 


N[jrj2A]  Q\,j2XP(X)]  RU1J2A,P(X)] 


(3-35a) 


Jl+J2 

N[j1,j2,X]=S-   (2j  +  l)I2j+1(4X) 

j=Jrj2 


(3-35b) 


Jl+J2  j+k 

I    (2j+l)U.  ,(4^)5 


P.7 


(3-35c) 


Q(jrj2,X,P(A.)]  =11  — l- 


J3J4  (2L +1)4  (2j4+l)4   2J3  2  2J42 


and 


N(jrj3,X]  N[j2j4,X]  NU3j4,X]  (3  ^ 


RUrj2A,P(?i)]  =  ?  1     1        1  4ui  (IP(X)) U5  (Ip(X)) 

J3J4  (2L +1)4  (2L+1)4   2J3  2  2J42 


N[jrj3,X]  N[j2,j4,X]  35 

The  B-curves  of  the  solutions  of  eq.  (2-35a)  with  different  truncations  in  j  are 
shown  in  Figure  10,  where  the  curves  again  converge  at  j=l  for  X  ranging 
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Table  3.  The  expectation  value  of  a  plaquette. 


analytical  calculation 

Monte  Carlo  simulation 

0.197(1) 

0.195871 

0.1964(08) 

0.262(1) 

0.259291 

0.2592(08) 

0.326(1) 

0.320871 

0.3194(08) 

-  0.393(2) 

0.384410 

0.3848(08) 

0.453(2) 

0.440293 

0.4404(08) 

0.516(4) 

0.498193 

0.4970(10) 

0.583(4) 

0.558514 

0.5572(12) 

0.646(3) 

0.614177 

0.6116(12) 

0.699(3) 

0.660372 

0.6632(12) 

0.763(4) 

0.715391 

0.7140(14) 

0.893(8) 

U.sl4o(14) 

0.996(5) 

0.907089 

0.9056(16) 

1.091(7) 

0.980752 

0.9810(16) 

1.186(7) 

1.050616 

1.0516(24) 

1.385(5) 

1.181302 

1.1830(10) 

1.554(5) 

1.272673 

1.2770(10) 
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Figure  7.  The  B-curves  of  the  expectation  value  of  a  plaquette  with  two 
Bianchi  identities  and  without  mean  plaquette  approximation. 


Figure  8.  The  B-curves  of  the  expectation  value  of  a  plaquette  with  12 
Bianchi  identities  and  without  mean  plaquette  approximation. 
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Figure  9.  The  convergent  B-curves  of  two  and  12  Bianchi  identities 

without  mean  plaquette  approximation,  and  that  of  12  Bianchi 
identities  with  mean  plaquette  approximation. 
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Figure  10.  The  B-curve  of  the  expectation  value  of  a  plaquette  with  12 
Bianchi  identities  and  with  mean  plaquette  approximation. 
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(a) 


(b) 


(c) 

Figure  11.  An  extended  planar  plaquette 

a)  Without  Bianchi  identities;  b)  With  four  Bianchi  identities; 
c)  With  20  Bianchi  identities. 
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from  0.05  to  1.1.  The  result  of  12  Bianchi  identities  with  the  mean  plaquette 
approximation  is  compared  in  Figure  9  with  the  results  of  two  and  12  Bianchi 
identities  without  the  mean  plaquette  approximation.  Since  the  largest 
difference  between  the  results  of  the  12  Bianchi  identities  with  and  without  the 
mean  plaquette  approximation  is  0.004  in  B-value  for  X  ranged  from  0.05  to 
1.1,  the  result  of  the  12  Bianchi  identities  with  mean  plaquette  approximation 
can  be  taken  as  the  convergent  result.  As  shown  in  Table  3,  this  result  agrees 
very  well  with  Monte  Carlo  calculation  results,  and  thus  the  effect  of 
nonabelianizability  of  Bianchi  identity  can  be  neglected. 

(b)  The  expectation  value  of  an  extended  plaquette:  To  calculate  the 
expectation  value  of  an  extended  planar  plaquette  which  contains  two  plaquettes 
(see  Figure  1  la),  the  proper  Bianchi  identities  must  be  built  up.  Four  Bianchi 
identities  can  be  used  to  cover  the  extended  plaquette  from  the  top  and  the  bottom 
(see  Figure  1  lb),  then  the  additional  16  Bianchi  identities  are  used  to  completely 
cover  the  surface  of  the  previous  four  Bianchi  identities  (see  Figure  1  lc).  Without 
mean  plaquette  approximation,  the  B-curves  of  the  four  and  the  20  Bianchi  identities 
with  different  truncations  in  j  are  shown  in  Figure  12  and  13,  the  curves  again 
converge  at  j=l  for  X  ranging  from  0.05  to  1.1.  To  combine  this  method  with  the 
mean  plaquette  approximation,  each  of  the  surface  plaquettes  of  the  Bianchi  identities 
must  be  replaced  by  its  expectation  value,  which  is  the  solution  of  eq.  (2-35a).  With 
20  Bianchi  identities,  the  expectation  value  is 
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i  z  l  i  cn  —L— )  Mij.ja]  Mu3,j  \]  N[uai 

1J1J2J3J4  1=1(2j.+l)4 

Fl  ^2       Z  Z  E  2  Z  (  n  • )  N[j  j  XJ  N[j  J  i]  NU  j  X] 

J1J2J3J4  1=1  (2j.+l)4 


N[y4,X]  S2tj1j2,j3,j4,X,P(X)]  T[j1,j2>X>P(X)]  T[j3J4,X,Pa)] 
N|jyj4,X]  S2U1,j2,j3,j4^,P(X)]  TOjjyX^CX)]  T[j3,j4,X,P(X)] 


TUrj3,?t,P(X)]  Tg2J£X£g:)] 
TUrj3,X,P(X)]  T[j2,j4A,P(X)] 


where 


4  U2k(2P(X)) 

SD1j2d3d4,x,P(X)]  =  z2:i2:  (  n  — s — — )NQ1,k1,X] 

klk2k3k4    1=1       (2k.+  l) 

N[jyk2,X]  N[j3,k3,X]  N[j4,k4,A.]  Nlkj.kyX] 
N[k3,k4,X]  N[krk3,X]  Nfkyk^X] 

and 

TTjlj^Pa)1=ZL(.n_i_^)N[jl,kl,M 

WzkjM  W&M  (3.36c) 
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The  B-curves  of  eq.  (2-36a)  are  also  convergent  at  j=l.  The  convergent  B- 
curve  of  eq.  (2-36a)  is  compared  in  figure  14  with  the  results  of  four  and  20 
Bianchi  identities  without  the  mean  plaquette  approximation.  Since  the 
difference  between  the  results  of  20  Bianchi  identities  with  and  without  the 
mean  plaquette  approximation  is  so  small  (0.01  for  the  largest  difference),  the 
result  of  20  Bianchi  identities  with  mean  plaquette  approximation  can  be  taken 
as  the  convergent  result. 


•1  .4  .7  1.0  1.3 

Figure  12.    The  B-curve  of  the  expectation  value  of  an  extended  planar 
plaquette  with  four  Bianchi  identities  and  without  mean 
mean  plaquette  approximation. 
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Figure  1 3.    The  B-curve  of  the  expectation  value  of  an  extended  planar 

plaquette  with  20  Bianchi  identities  and  without  mean  plaquette 
approximation. 

3.2.2  Discussions  and  Comments 

The  analytical  calculation  method  presented  above  is  more  than  simply  an 
alternative  calculational  method.  Rather  it  reveals  insights  into  the  underlying 
physics  and  physical  structure  of  the  lattice  system.  For  example,  from  the 
change  in  B-curve  that  occurs  as  a  function  of  the  number  of  the  Bianchi 
identities,  the  effect  of  finite  lattice  size  can  be  understood  explicidy  and  thus  the 
lattice  size  needed  to  avoid  the  finite-size  effect  can  be  estimated.  The 
B-curve  also  contains  some  important  information  about  the  lattice  system  and  the 
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Figure  14.    The  convergent  B-curves  of  four  and  20  Bianchi  identities  without 
mean  plaquette  approximation,  and  that  of  20  Bianchi  identities 
with  mean  plaquette  approximation. 

trial  wavefunction  proposed  in  the  calculation.  A  typical  B-curve  of  an  arbitrary 
physical  observable  is  shown  in  Figure  15.  The  curve  can  be  divided  into  three 
regions:  One  is  for  g  ^  gj»  B-value  approaching  zero,  and  with  the  system  in 
the  strong  coupling  region;  another  is  for  %2  <  8  <  Sj>  B  rising  rapidly,  and 
with  the  system  in  the  intermediate  coupling  region;  and  the  third  is  for  g  <  g^, 

B-value  remaining  a  constant,  and  with  the  system  in  the  weak  coupling  region. 
As  seen  in  Figure  9,  this  is  not  avoidable  for  the  decline  in  the  B-curve  for 


Figure  15.  A  typical  B-curve  as  a  function  of  coupling  constant. 


X  >  0.9  ( — !—  >  0.38  )  even  if  the  number  of  the  Bianchi  identities  is  increased. 
This  situation  indicates  that  the  trial  wavefunction  proposed  is  not  good  enough 
to  incorporate  sufficient  correlations  in  the  weak  coupling  region. 
Unfortunately,  the  analytical  method  presented  above  is  effective  only  when  the 
trial  wavefunction  is  localized.  For  example,  consider  a  two-parameter  trial 
wavefunction: 


¥  =  exp{X  XtrU+A..  ItrUJ 

1  P     p     2  s  P 
V  P  (3-37a) 


where 


U_  =  the  extended  plaquette  shown  in  Figure  16(b) 
p  (3-37b) 
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To  calculate  the  expectation  value  of  plaquette,  we  first  evaluate  the  function 


=  ±-  In  { J  [  PI  dU,  ]  exp  [  2X.ZtrU_  +  2X,  ZtrU_  ]  } 
IN.,  k      *  p     v         p  p 


=  ^-{ln(J[ndUk]exp[2X1ZtrUp]) 


J  [  n  dU.  ]  exp  [  2X,  I  trU  +         trU_  ] 

+  ln(      k    fk  1P      P       ^      P  )} 

J  [  n  dU,  ]  exp  [  2^  E  trU_  ] 
k     k  2p  p 

-j-{ln(J[ndU.  ]exp[2X1ZtrU_])  +  ln<exp(2X2ZtrU_)>} 
JN  k      k  xpP  pP° 


-  ^  ( in  ( J[  n  dUk  ]  exp  [  2X,  S  trUp  ]  >  +2  <  (  ?  TU,  )"  >  ) 

where  Np  is  the  total  number  of  plaquettes;  <  V  ^  is  the  expectation  value  of 

physical  observables  V,  obtained  by  using  the  wavefunction  of  eq.  (3-37a)  with 
^2  =0,  and  <  Vn  >  stands  for  the  connected  part  with 


<v>c  =  <v>0 


<Vn>  =  <Vn^-<V^  forn>l 


(3-39b) 
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The  expectation  value  of  a  plaquette  then  is 


2  dX 


1 


n 


=  <  trLL  >  + 


<„5 


l  n! 


<(ZtrUJn>} 

n         D  v 


P  P 


(3-40) 


the  second  term  on  the  right  hand  side  is  not  easily  manageable  for  n  >  1  using 
the  analytical  method.  However,  the  product  of  single  plaquette  wavefunction, 
eq.  (3-28),  is  exact  in  the  strong  coupling  region  and  the  scaling  region  for 
SU(2)  gauge  theory  is  in  the  range  of  0.25  <       <  0.3  (0.3  <  X  <  0.7),  it  is 


thus  reasonable  to  use  eq.  (3-28)  as  the  trial  ground  state  to  carry  out  the 
calculation  down  to  the  scaling  region.  Whether  or  not  the  continuum  physics 
can  be  extracted  from  the  calculation,  it  is  necessary  to  check  if  the  result  obeys 
the  scaling  behavior  as  discussed  in  Chapter  2.3. 

To  conclude  this  section,  final  points  should  be  noted: 

(a)  From  the  two  examples  of  calculations  worked  out  here,  a  general  rule 
for  building-up  a  sufficient  number  of  Bianchi  identities  can  be  developed.  For 
a  physical  observable  which  is  a  function  of  plaquettes,  enough  Bianchi 
identities  first  must  be  used  to  cover  completely  the  considered  plaquettes  from 
top  and  bottom.  Then  additional  Bianchi  identities  are  used  to  cover  completely 
the  surface  plaquettes  of  the  previous  Bianchi  identities.  To  perform  the 
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calculation,  the  surface  plaquettes  of  the  whole  Bianchi  identities  established  are 
replaced  by  the  mean  plaquette  and  all  of  the  Bianchi  identities  are  truncated  at 
j=l.  This  general  rule  will  be  applied  in  carrying  out  the  calculations  that 
follow. 

(b)  Considering  both  planar-  and  a  nonplanar-extended  plaquettes  (see 
Figure  16),  one  Bianchi  identity  is  needed  to  correlate  two  plaquettes  in  a 
nonplanar-extended  plaquette  and  two  Bianchi  identities  are  needed  for  the  case 
of  a  planar-extended  plaquette.  By  using  four  Bianchi  identities,  the  B-curves 
of  the  expectation  values  of  these  two  quantities  are  shown  in  Figure  17.  From 
the  B-curves  it  is  possible  to  see  that  the  effect  of  the  Bianchi  identities  on  the 
nonplanar  case  is  larger  than  the  effect  on  the  planar  case.  In  a  similar  manner, 
one  can  determine  the  most  probable  configuration  for  the  collective  mode 
which  is  crucial  in  the  study  of  the  glueball  mass. 


(b) 


Figure  16.  The  extended  plaquettes 

a)  Planar;  b)  Non-planar. 
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Figure  17.  The  B-curves  of  the  expectation  values  of  a  planar-  and  a 
nonplanar-extended  plaquette  with  four  Bianchi  identities. 


CHAPTER  4 

VARIATIONAL  VACUUM  AND  STATIC  QUARK  POTENTIAL 
FOR  SU(2)  HAMILTONIAN  LATTICE  GAUGE  THEORY 


The  calculation  of  the  vacuum  energy  is  most  natural  in  the  Hamiltonian 
formalism  since  the  vacuum  energy  is  the  ground  state  expectation  value  of  the 
lattice  Hamiltonian 

<y0IHiyQ> 
0  <VQ\VQ> 

The  vacuum  energy  is  not  a  physical  quantity.  However,  the  difference 
between  the  vacuum  energy  density  in  the  strong-  and  the  weak-coupling  limits 
may  have  some  physical  implications  since  this  is  just  the  "Bag"  constant  in  the 
M.I.T.  Bag  model  (*)  of  hadrons.  The  crucial  step  in  any  variational 
calculation  is  to  propose  a  reasonable  trial  ground  state  wavefunction  that  is 
capable  of  describing  the  vacuum  properties.  This  will  be  the  subject  of 
Section  1. 

The  static  quark  potential  is  a  physical  quantity  which  is  related  to  the 
vacuum  structure.  Introducing  a  quark-antiquark  pair  into  the  vacuum, 
accomplished  in  Chapter  2,  permits  study  of  the  static  quark  potential  which  is 
defined  as  the  energy  difference  between  the  vacuums  with  and  without  the 
fermion  pair. 
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4.1.  Ground  State  Energy  per  Plaquette 
The  Kogut-Susskind  Hamiltonian  is 


aH  =  -^X2E?+44l[2-trU  ] 
4%  k  a  k       p  P 


(4- la) 


where  £  =  — ^  ,  and  the  Ea  arc  the  color-electric  field  operators  which  satisfy 
2g  K 


the  commutation  relation 


^  V=\.c-fUk 


where  X    are  the  Pauli  matrices.  The  trial  ground  state  wavefunction  used  is 


¥0  =  exp{X2trUJ 

p  (4-2) 


where  X  is  a  variational  parameter.  As  discussed  in  Chapter  3.2.2,  this  trial 
wavefunction  is  exact  in  the  strong-coupling  limit,  and  it  is  a  reasonable 
approximation  in  the  intermediate  region.  Using  integration  by  parts,  it  can  be 
shown(37)  that 


<  To  1 1 lK K '  * 0  >4 <W0  1  g  '  K  Ek  • *0  11 '  *0  >  (4.3a) 
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Thus,  the  kinetic  energy  is 


<  *P,J  XI  E?  E?  I  > 


0  ka   k  k     0    _3.N  pa) 

0     0  (4-3b) 


where  Np  is  the  total  number  of  plaquettes  and 


<  *¥n  1 1  trUn  I  Wn  > 


P  0     0  (4-3c) 


is  the  expectation  value  of  a  plaquette.  The  ground  state  energy  per  plaquette  is 


1    <Y  laHI*P  >  . 

X>P     <Y0IY0>         8^  (4.43) 


and  the  variational  parameter  ^  is  deteimined  by  minimizing  a£Q  with  respect  to 
X,  which  gives 


128    4p(*)  (4-4b) 


where  p'(X)  is  the  derivative  of  P(A.)  with  respect  to  X.  Therefore,  to  find  the 
ground  state  energy  per  plaquette,  the  quantity  that  must  be  calculated  is  the 
expectation  value  of  a  plaquette.  Using  the  result  obtained  in  Chapter  3.2. 1, 
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where  the  expectation  value  of  a  plaquette  is  calculated  by  using  12  Bianchi 
identities  with  mean  plaquette  approximation,  the  resultant  ground  state  energy 
per  plaquette  is  shown  in  Figure  18. 


.2  .3  .4 

Figure  18.  The  ground  state  energy  per  plaquette  as  a  function  of  coupling 
constant. 

The  exact  ground  state  energy  per  plaquette  also  can  be  obtained  through 
numerical  computation.  The  utility  of  applying  stochastic  techniques  developed 
for  quantum  many-body  systems  to  the  study  of  Hamiltonian  lattice  gauge 
theory  has  been  demonstrated.  Chin  et  al.(38)  applied  the  guided-random  walk 
algorithm  to  solving  the  ground  state  energy  of  the  SU(2)  lattice  gauge  theory. 
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In  Table  4,  our  result  for  the  ground  state  energy  per  plaquette  is  compared  with 
the  Monte  Carlo  calculations.  The  second  column  in  Table  4  displays  Chin  et. 
al.'s  Monte  Carlo  result(38)  using  the  variational  ground  state  eq.  (4-2).  As  can 
be  seen  our  analytical  calculation  agrees  very  well  with  the  Monte  Carlo 
calculation.  The  third  column  in  Table  4  provides  the  exact  results  obtained  by 
Chin  et.  al.,(38)  the  difference  between  the  variation  results  and  the  exact  ones 

appear  only  after  the  second  decimal  for  the  scaling  region,  which  is  about  the 

range  of  0.25  <-hr<  0.3.  It  therefore  can  be  seen  that  the  vacuum  fluctuation 
2g2 

is  quite  localized  and  that  the  magnetic  fluctuations  are  unlikely  to  correlate 
plaquettes  separated  by  more  than  one  unit  of  lattice  spacing  for  the  intermediate 
coupling  region.  If  the  magnetic  fluctuations  of  two  units  of  lattice  spacing  can 
be  included  in  the  trial  ground  state,  the  difference  between  the  variational 
results  and  the  exact  ones  will  be  reduced  for  the  intermediate  coupling  region. 
However,  because  of  the  complexity  involved  in  eq.  (3-40),  eq.  (4-2)  will  be 
used  as  the  trial  ground  state  to  perform  analytical  calculations  for  the  static 
quark  potential  and  the  glueball  mass  of  SU(2)  Hamiltonian  lattice  gauge 

theory. 

4.2.  Static  Quark  Potential 

To  give  a  clear  global  physical  picture  to  the  lattice  calculation  of  the  static 
quark  potential,  a  simple  model  is  presented  in  the  first  part  of  the  section.  The 
model  is  that  of  a  rigid  vacuum  with  fluctuating  string,  which  originally  was 
constructed  by  Liischer.(39)  In  the  second  part,  the  analytical  method 
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Table  4.  The  variational  and  exact  ground  state  energies  per  plaquette. 


analytical  and  variational 

Monte  Carlo  and  variational 

exact 

0.136(9) 

0.28510093 

0.28510(001) 

0.28510(006) 

0.158(1) 

0.37362619 

0.37365(001) 

0.37356(011) 

0.176(8) 

0.45902315 

0.45903(002) 

0.45872(006) 

0.193(6) 

0.54136152 

0.54136(002) 

0.54050(009) 

0.209(2) 

0.62072482 

0.62073(004) 

0.61943(011) 

0.223(6) 

0.69720007 

0.69727(006) 

0.69527(013) 

0.237(2) 

0.77085766 

0.77079(005) 

0.76742(022) 

0.25 

0.84176066 

0.84171(007) 

0.83581(028) 

0.262(2) 

0.90997328 

0.90989(008) 

0.90231(032) 

0.273(9) 

0.97556613 

0.97560(010) 

0.96682(042) 

0.295(8) 

1.09917318 

1.09913(020) 

1.08357(064) 

0.316(2) 

1.21321877 

1.21305(023) 

1.18809(058) 

0.335(4) 

1.31859769 

1.31875(020) 

1.28321(083) 

0.353(6) 

1.41637802 

1.41688(028) 

1.37865(102) 

0.395(3) 

1.63488530 

1.63615(033) 

1.57650(079) 

0.433(0) 

1.82695077 

1.82586(070) 

1.76400(112) 
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developed  is  used  to  calculate  the  static  quark  potential.  Using  the  static  quark 
potential  as  a  function  of  the  separation  distance  between  quark  and  antiquark, 
the  critical  coupling  strength  for  the  restoration  of  rotational  symmetry  is  located 
and  the  string  tension  is  extracted  by  fitting  the  potential  to  a  simple  polynomial 
of  the  separation  distance.  Finally,  the  scaling  behavior  of  string  tension  is 
examined  by  comparing  this  behavior  with  the  expected  asymptotic  scaling  line. 

4.2.1.  Simple  Model  for  Static  Quark  Potential:  Rigid  Vacuum  with  a 
Fluctuating  String 

From  the  simple  calculation  using  strong-coupling  expansion  in  Chapter 

3.1,  one  can  see  that  the  chromoelectric  flux  tube  is  formed  in  the  strong 

coupling  limit.  Because  of  the  potential  energy  term,  ^2  p  ^ "  ^P^'  ^e 
flux  tube  starts  to  fluctuate  when  the  coupling  strength  decreases.  On  the  other 
hand,  it  was  demonstrated  that  the  magnetic  fluctuations  in  the  vacuum  are  quite 
localized.  Combining  these  two  facts,  we  assume  that  the  vacuum  is  rigid,  and 
in  the  presence  of  a  static  quark-antiquark  pair  at  a  finite  distance,  L,  the  glue 
field  can  be  described  by  a  fluctuating  thin  tube  containing  the  chromoelectric 
flux  which  flows  from  the  quark  to  the  antiquark.  Since  interest  here  is  only  in 
the  low-momentum  mode,  the  degrees  of  freedom  inside  the  flux  tube  are  not 
disturbed.  Thus,  the  dynamics  of  a  relativistic  thin  string  can  be  considered, 
and  the  result  is  independent  of  the  gauge  group. 

Supposing  that  the  quarks  are  located  at  x  =  ( 0, 0,  0)  and  x  =  ( L,  0,  0)  ^ 
the  string  then  can  be  described  by  a  two-component  vector  field 
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lu^x1),  OSx^Uwith  4(x°,0)  =  ^(x°,L) 


as  shown  in  Figure  19.  To  study  the  dynamics  of  a  thin  string,  the  effective 


Figure  19.  The  ^-field.defined  in  the  dynamics  of  a  thin  string. 

action  for  the  ^-field.  must  be  known.  As  discussed  in  the  beginning  of 
Chapter  2,  the  symmetry  manifested  by  the  system  is  the  most  plausible 
guideline  for  constructing  the  action.  Therefore,  the  action  can  be  expected  to 
have  the  following  general  properties: 

(a)  the  effective  action  must  be  local,  i.e. 


A 


q 


1 


(4-5) 


where  the  Lagrangian  density,  L(x),  depends  on  the  ^-field.  and  its  derivative. 
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(b)  L(x)  should  be  invariant  under  the  following  symmetry  operations: 

(i)  Pointcare  transformation  in  the  ( x* )  plane  ; 

— ♦ 

(ii)  0(2)  rotations  and  translation  of  the  £  field. 

-•  -* 

The  property  (b-ii)  precludes  mass  term  £(x)  •  ^(x)  in  the  Lagrangian  density. 
The  string  thus  is  massless  and  L(x)  must  be  made  up  of  derivatives 


The  effective  action,  eq.  (4-6),  is  not  renormalizable  and  the  parameters  b  and  c 
are  dimensional.  However,  only  the  low-momentum  mode,  where  the 
wavelength  is  much  longer  than  the  intrinsic  width  of  the  string,  is  of  interest 
here.  For  the  high-momentum  mode,  the  internal  degrees  of  freedom  of  the 
tube  are  excited,  and  the  physical  picture  expressed  by  the  fluctuating  string 
breaks  down.  Therefore,  those  terms  with  dimensional  parameters  are 
irrelevant  to  low-energy  physics  and  it  is  sufficient  to  take 


3^,3^$,  -,  (n,  -u  =  0,  1). 


So  as  a  result 


+  c(3^/?)2  +•••} 


(4-6) 
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s   =j  dxOdx^aJ./u 

ett     0<xx<L  K  (4-7) 

as  the  effective  action.  There  are  two  important  features  of  this  action: 

(a)  The  roughening  effect  from  the  fluctuation  in  the  string:  Because  the 

K  ~*        1  1 

£  -  field  is  massless,  the  variance  in  ^(xQ,x  )  for  any  x    between  0  and  L  is 


<i2(xo,xl)>.JA.ln(L) 


where  a  is  the  intrinsic  width  of  the  tube.  Since  the  variance  is  proportional  to 
'  an  infinitely  long  localized  flux  tube  cannot  exist  in  a  continuum 

Yang-Mills  theory.  In  lattice  gauge  theories,  the  situation  is  quite  different.  In 
this  case,  the  translational  symmetry  is  only  discrete.  Consequently,  there  is  a 
mass  gap  which  is  proportional  to  g2,  separating  the  excited  mode  and  the 
ground  state.  For  a  coupling  strength  just  below  the  critical  value  g£  ,  the  mass 
gap  becomes  easy  to  overcome,  the  flux  tube  appears  rough,  and  the  bulk 
rotational  symmetry  is  restored. 

(b)  The  correction  to  linear  confinement:  Because  of  the  zero-mass 
fluctuation  of  the  string,  there  are  power  corrections  to  the  linear  quark  potential 


V(L)  =  aL  +  ^  +  -|-  +  0(L3) 

L  (4-9) 
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where  a  is  the  string  tension,  a  is  a  dimensionless  parameter,  and  p  is  a 
parameter  with  dimension  of  length.  The  ^  term  is  due  to  a  one-dimensional 

Casimir  effect(22)  and  the  coefficient  is  universal 
a  =  -JL»-0.26 


where  the  minus  sign  implies  the  force  is  attractive.  In  the  lattice  gauge  theories, 
because  of  the  existence  of  the  mass  gap,  eq.  (4-9)  is  modified  as 


V(L)  =  gL  +  £  +  £-  +  0(  L"2) 

a     L  (4-10) 


where  the  constant  5  is  dimensionless  and  corresponds  to  the  mass  gap  between 
the  ground  state  and  first  excited  state  of  the  flux  tube.  The  string  tension  a  has 
dimension  of  mass  squared  and  is  a  function  of  coupling  constant.  In  the  next 
section,  a  will  be  determined  by  calculating  V(L)  at  various  values  of  L  and 
coupling  constant  Sq  . 

4.2.2.  Static  Quark  Potential  for  SIK2)  Hamiltonian  Lattice  Gauge  Theory 

In  this  section  where  the  static  quark  potential,*  is  calculated,  it  is 
necessary  to  impose  a  quark-antiquark  pair  into  the  vacuum.  Consider  a 


*For  the  calculation  of  the  static  quark  potential  in  two  spatial  dimensions,  see 
ref.  (40). 
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qq  state,  where  q  and  q  denote  a  heavy  quark-antiquark  pair,  respectively:  Let 
q  and  q  be  created  by  the  operators  ^(n)  and  4/+(n).  The  generator  of  a  color 
rotation,  which  was  defined  in  eq.(2-27),  now  is  generalized  to 

I  E*  (n)  -»  X      (n)  +  XP+  (n)^¥(n) 

k    K  k  1  (4-10) 

at  site  n.  To  describe  the  heavy  qq  state  one  should  choose  an  operator  which 
is  invariant  under  the  transformation  (4-10).  The  operator 


¥+(n){  n    U.  }¥(n  +  R) 

kepath   k  (4Ua) 


where  the  path  extends  from  n  to  n  +  R  but  otherwise  is  arbitrary,  satisfies  this 
requirement.  Combining  the  operator,  eq.  (4- 11  a),  with  the  vacuum  obtained  in 
Section  4.1  of  this  chapter,  the  qq  state  in  the  gluon  field  background  is 
described  by 


(4-1  lb) 


where  I  ¥Q  >  is  the  variational  vacuum  defined  in  eq.  (4-2).  Note  that  the 
vacuum  of  eq.  (4-1  la)  is  rigid  and  independent  of  the  coupling  strength.  In  the 
strong-coupling  region,  there  is  no  significant  difference  between  eqs.  (4-1  la) 
and  (4-1  lb)  because  the  variational  vacuum  is  almost  rigid,  i.e.  the  variational 
parameter  X  of  eq.  (4-2)  is  small  in  the  strong  coupling  region.  But  the 
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difference  is  significant  in  the  intermediate  region  where  localized  magnetic 
fluctuations  exist  in  the  variational  vacuum. 

For  a  quark  and  antiquark  pair,  the  quark  is  placed  at  the  origin  and  the 
antiquark  is  at  lattice  site  (nx,nv,nz).  In  general,  there  are  infinite  ways  to 
choose  a  path  to  go  from  the  origin  to  the  lattice  site  (nx,nv,nz).  The  dimension 
of  the  Hilbert  space  spent  by  the  states  of  eq.  (4-1  lb)  thus  is  infinite.  To  pick 
up  a  subset  of  the  Hilbert  space,  the  following  two  important  simplifications  are 
made: 

(a)  The  paths  chosen  in  eq.  (4-1  lb)  are  those  which  possess  the  shortest 
distance  from  the  origin  to  the  lattice  site  (nx,ny,nz).  The  states  thus  defined  in 
eq.  (4-1  lb)  give  the  lowest  energy  in  the  leading  order  of  strong-coupling 
expansion. 

(b)  To  further  reduce  the  number  of  states,  the  antiquark  is  always  placed 
in  the  x-y  plane. 

After  these  two  simplifications,  each  path  in  eq.  (4-1  lb)  contains  the  (nx+ 

(n  +  n  )! 

ny)  links,  and  the  dimension  of  this  subset  is  —  £—  for  a  particular  set  of 

nx!  ny! 

nx  and  ny. 

After  choosing  the  subspace,  the  next  step  is  to  compute  the  matrix 
elements  of  the  Hamiltonian  with  respect  to  the  states  in  the  subset  Let  the  ith 
state  be 


li>  =  y+(„x,ny)F.>P(0,0)l¥o>  ^ 

and 

F.  =    n  U, 

i  kepathi    *  (4- 12b) 


where  is  given  in  eq.(2-3).  The  contribution  from  the  magnetic  potential  to 
the  matrix  element  of  the  Hamiltonian  between  states  I  i  >  and  I  j  >  is 


<ilZ(2-trU  )lj> 

 2  1  =  2Nn<tr(F^F.)>-l^<tr(F^F.)> 


-  < I  trU  >< tr  ( F* F. )  > 
p     P  i  j' 


<E(2-trU  )><tr(F*F  )  >-l  JL<tr(F+F )> 
P  P  1  J       2  dX  1  J  (4-13) 


The  left  hand  side  of  eq.  (4-13)  is  divided  by  <Vq  I  Vq>  for  convenience  and 
the  expectation  values  on  the  right  hand  side  are  evaluated  with  respect  to  the 
variational  vacuum  I      >  in  eq.  (4-2).  The  contribution  from  the  chromo- 
electric  term  in  the  Hamiltonian  is 
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<ilILE*E*lj> 

 ka   k   =  1\  { I^<tr(F^F.)>  +  <2:trUn><tr(F^F.)> 

<v0 1  vQ>       2     2  ax       1  J        P    P         1  J 


+  |(n  +n  +C  (n  n  ))<tr(F+F  )> 

8    x    y    ij  x  y         i  j  (4-14) 


where  C.j(nx,riy)  is  the  number  of  common  links  between  Fi  and  Fj.  The 
derivation  of  eq.  (4-14)  is  shown  in  Appendix  A.  Combining  eqs.  (4-13)  and 
(4-14)  yields 


<ilaHlj>=  +R 
<V0%>      P   8$  1  J 


+  (Jk_2^)A-<tr(F*F.)> 
16^        BX         i  J 


+  3%[nx+ny+cij(nx'V]<tr(Fi"Fj)> 


(4-15) 


where  the  first  term  on  the  right  hand  side  of  eq.  (4-15)  is  simply  the 
contribution  from  the  vacuum  energy  defined  in  eq.  (4-4a).  Therefore,  the 
matrix  element  of  the  static  quark  potential  between  states  I  i  >  and  I  j  >  is 

Vij(nx,ny)  =  +  (^-20^<tr(FirFj)> 


+  l2V[nx  +  ny  +  Cij(nx'V]<tr(Fi"Fi):> 

^  J  J  (4-16) 
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Since  the  states  which  belong  to  the  chosen  subset  are  not  orthogonal  to 
each  other 


M..5   <1|J>  =<tr(F^F.)> 

y  <v0i¥0>        1  j 


(4-17) 


As  a  result,  it  is  necessary  to  diagonalize  the  matrix  M  to  form  a  set  of 
orthogonal  states  and  then  to  normalize  each  state.  This  is  just  the  Schmidt 
orthogonalization  process  and  the  process  yields  a  set  of  orthogonal  states 
{  I  k  >  }  with 


k>  =  (SD  \.\i> 

ki  (4-18) 


where  S  is  the  transformation  which  diagonalizes  the  matrix  M, 


S+MS  =  D 


(4-19) 


l 

D  is  a  diagonal  matrix,  and  D  2  also  is  a  diagonal  matrix  with  diagonal 
elements  equal  to  the  inversed  square  root  of  the  diagonal  elements  of  D.  In 
terms  of  orthogonormal  states  {  I  k  >  } ,  the  matrix  element  of  the  static  quark 
potential  between  states  I  k  >  and  I  k'  >  is 


vkk,  =  (D  V)    V  (SD  ~2) 

a   y  JK  (4-20) 
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By  diagonalizing  the  matrix  V,  the  eigenvalues  are  obtained  and  the  lowest 
eigenvalue  gives  an  upper  boundary  on  the  ground  state  energy  of  the  heavy  qq 
state. 

N„ 


i— < — ■ — <- 


>    ■  > 


-►  N, 


Figure  20.  An  example  of  tr(Fi+Fj)  for  nx  =  2,  ny  =  1,  and  Cij  =  1. 


The  quantities  that  must  be  calculated  to  obtain  the  matrix  elements  of  M 
and  V  are  <  tr(Fi+Fj)  >  for  all  i  and  j.  An  example  of  tr(Fi+Fj)  for  nx  =  2, 
ny  =  1  is  shown  in  Figure  20.  Since  the  quark  is  placed  at  the  origin  and  the 
antiquark  is  always  in  the  x-y  plane,  tr(Fi+Fj)  is  the  trace  of  a  planar  extended 
plaquette  for  any  i  and  j.  Because  of  the  simplicity  of  the  variational  vacuum  in 
eq.  (4-2),  it  is  always  possible  to  use  the  Gross- Witten  trick  (21)to  convert  the 
expectation  value  of  the  trace  of  an  extended  plaquette  to  the  expectation  value 
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of  the  product  of  the  traces  of  plaquettes  contained  in  the  extended  plaquette. 
To  explain  the  Gross-Witten  trick,  consider  the  expectation  value  of  <|>(P) 
where 

<|)(P)  =  tr  (the  extended  plaquette  shown  in  Figure  16(a)) 

=tr(uPlV 

Pl  P2  (4.2l) 
Since  the  variational  vacuum,  eq.  (4-2),  is  a  function  of  trUp  only,  it  is 
invariant  under  the  transformation 


u_-»vu_v+ 

P         P  (4-22) 


where  V  is  an  element  of  SU(2)  gauge  group.  The  expectation  value  of  <|>(P) 
then  is 


J  [  n  dU  ]  [  tr(UD  UD  )  ]  exp(  2X  Z  trU_  ) 

<<D(P)>=    k   rk  P-±-h  P  P 

J  [  n  dU,  ]  exp(  2X  £  trU  ) 
k     K  p  P 


J  [  n  dUk  ]  [  dV  ]  [  tr(VUp  V+Up  )  ]  exp  {  2\  Z  trUp  } 

JtridU,  ]exp{  2XZtrUn  } 

k     k  p     p  (4_23) 
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For  SU(N)  group 


+  =J_ 

jk    N  "ik  1j  (4_24) 


JdVVu\ 


It  then  follows  from  eqs.  (4-23)  and  (4-24)  that 


<<Kp)>=y<trU    •trU  > 

1        v\        v2  (4-25) 


The  process  of  eq.  (4-23)  can  be  repeated  to  obtain  the  following  general  result: 
For  an  extended  plaquette  f(p)  containing  m  plaquettes,  U  ,  U  ,  •  •  • ,  U  , 
by  applying  the  Gross- Witten  trick(21)  repeatedly  it  can  be  shown  that 


1  m 
<trf(p)>  =  -±—  <n  trU  > 

2         i=l       1  (4-26) 


Therefore,  in  the  lattice  calculation  of  static  quark  potential,  the  quantities  that 
must  be  calculated  are  the  expectation  values  of  various  products  of  the  traces 
of  plaquettes.  In  summary,  the  variational  parameter  X  is  determined 
by  minimizing  the  ground  state  energy,  eq.(4-4a),  then  the  potential, 
eq.(4-16),  is  evaluated  by  using  the  Gross-Witten  trick(21)  and  the  technique 

developed  in  Chapter  3.2.  The  results  of  static  quark  potentials  for 
(nx=l,ny=3),  (nx=2,ny=2),  (nx=2,ny=3),  and  (nx=3,ny=3)  as  functions  of 
coupling  constant  are  shown  in  Figure  21. From  the  results  of  Figure  21,  one 
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can  locate  the  roughening  transition  point  and  extract  the  string  tension  for 
SU(2)  Hamiltonian  lattice  gauge  theory: 

(a)  The  roughening  transition  point:  In  Section  2. 1  of  this  chapter,  it  was 
pointed  out  that  in  lattice  gauge  theories  a  mass  gap  exists  for  the  fluctuating 
flux  tube  and  that  mass  gap  easily  can  be  overcome  when  the  coupling  strength 
decreases.  Once  the  mass  gap  can  be  overcome,  the  flux  tube  starts  to  fluctuate 
and  the  rotational  symmetry  effectively  is  restored.  In  Figure  22,  the 
equipotential  surfaces  of  the  static  quark  potentials,  shown  in  figure  21,  are 
drawn  for  various  ^  values.  The  rotational  symmetry  is  restored  at  JL  =  0.26 


Figure  21.  The  results  of  static  potentials  for 

(nx=1'ny=3>'  (nx=2,ny=2),  (nx=2,ny=3),  and  (nx=3,ny=3). 
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which  is  the  critical  coupling  strength  for  the  roughening  transition.  It  remains 
rotationally  symmetrical  from  £  =  0.26  up  to  £  =  0.32  After  that,  the  perfect 
spheres  start  to  be  distorted  slowly.  Two  statements  should  be  made  about  the 
critical  coupling  strength  ^c  =  0.26  : 

(i)  The  result  of  the  calculation  applied  here  is  very  close  to  the  value  obtained 
by  Monte  Carlo  calculation  in  Lagrangian  formalism  where  £c  =  0.22,  but  quite 
different  from  the  strong-coupling  expansion  result  (41)  in  the  Hamiltonian 

formalism  where  P  =  o.5-  As  discussed  in  Chapter  2.2.2,  for  the  same 
c 

physical  phenomena  g    <  g     However  the  ratio  of  Au  and  AT  is  very  close 
rl       L,  ii  L, 

to  1,  £    in  the  Hamiltonian  formalism  is  expected  to  be  slightly  higher  than 
c 

F    in  the  Lagrangian  formalism.  On  the  other  hand,  for  the  Hamiltonian 

c 

calculation  the  vacuum  in  the  strong-coupling  expansion  is  assumed  to  be  rigid, 
and  the  variational  vacuum  allows  magnetic  fluctuation.  Therefore,  £ 

c 

calculated  by  strong-coupling  expansion  is  expected  to  be  higher  than  the  result 
of  the  variational  calculation. 

(ii)  Recent  Monte  Carlo  calculations(31)  of  the  glueball  mass  in  the 
Hamiltonian  formalism  show  that  the  strong-coupling  expansion  cannot  give  the 
correct  result  when  £  >  £0  «■  0.3  for  SU(3).  The  £0  should  be  slightly  smaller 

than  0.3  for  SU(2),  because  the  correlation  effect  of  SU(2)  is  stronger  than  that 
of  SU(3).  Since  the  glueball  mass  is  inversely  proportional  to  the  fluctuated 
width  of  the  flux  tube  and  the  reason  for  the  strong-coupling  expansion 
breaking  down  at  the  rough  transition  point  in  the  calculation  of  string  tension 
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was  already  explained  in  Chapter  3. 1,  we  therefore  conclude  ^  =  which  is 
consistent  with  the  result  we  obtained. 

(b)  The  string  tension:  To  find  the  string  tension,  the  static  quark  potential 
is  fitted  to  a  simple  form 


V(L)a  =  (aa2)^  +  5  +  -^- 

a         (L/a)  (4-27) 


where  a  is  the  lattice  spacing,  a  is  the  string  tension,  the  constant  5  is  due  to  the 
appearance  of  mass  gaps  in  the  lattice  version  of  the  string,  and  the  Coulomb 
potential  is  due  to  a  one-dimensional  Casimir  effect  Using  the  least-squares 
fitting,  the  string  tension  is  extracted  from  the  static  quark  potentials  between  a 
quark  located  at  the  origin  and  an  antiquark  at  (nx  =  1,  nv  =  3), 
(nx  =  2,  ny  =  2),  (nx  =  2,  nv  =  3),  (nx  =  3,  ny  =  3),  respectively,  with  the 
results  shown  in  Figure  23.  The  data  points  of  the  string  tension  are  shown 
only  up  to  £  =  0.32  because  the  variational  vacuum  is  not  accurate  after  that 
point.  The  scaling  behavior  exhibited  by  the  string  tension  agrees  very 
well  with  the  expected  asymptotic  scaling  line  which  yields 


ATT  =  0.009  Jo 


(4-28) 
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\  =  0.32  4  =  0.4 


Figure  22.  The  equipotential  surfaces  for  various  values  of  %. 
The  small  dots  denote  lattice  sites. 
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.11  .15  .19     „       .23  .27  .31 


Figure  23.  The  string  tension  as  a  function  of  the  coupling  constant. 
The  solid  line  is  the  expected  asymptotic  scaling  line. 


If 


A.  =(0.01 19  ±0.0015) /a 
L  (4-29) 


which  is  the  best  estimation  from  the  Lagrangian  calculation.is  adopted,^)  the 
ratio  yielded  is 


104 


H=  0.76  ±0.11 


\  (4-30) 


which  is  quite  close  to  the  theoretical  prediction 
AH 

-S-  =  0.84 

^  (4-31) 


obtained  by  Hasenfratz  and  Hasenfratz.(33) 


CHAPTER  V 
GLUEBALL  MASS  FOR  SU(2)  HAMHLTONIAN 
LATTICE  GAUGE  THEORY 


A  color  singlet  bound  state  of  color  gauge  bosons  is  called  a  glueball  state.  In 
terms  of  string  physics  which  were  discussed  in  Chapter  4.2.1,  glueball  mass  is 
inversely  proportional  to  the  fluctuating  width  of  the  flux  tube,  Thus,  the  glueball 
state  is  an  excited  collective  mode  of  the  system.  Glueball  masses  with  different 
jPC  assignments  represents  the  energy  differences  between  these  low-lying 
excited  states  and  the  ground  state.  This  chapter  is  devoted  to  calculations  and 
discussions  of  0++  glueball  mass  for  SU(2)  Hamiltonian  lattice  gauge  theory. 

To  calculate  the  glueball  mass,  the  variational  vacuum—obtained  in 
Chapter  4.1 -is  used  as  the  ground  state,  and  the  trial  lowest  excited  state  is 
constructed  to  include  correlations  of  up  to  two  adjacent  plaquettes.  Calculations 
and  results  are  presented  in  Section  5.1.  Section  5.2  is  reserved  for  discussion  and 
comments  regarding  the  results  obtained  in  Section  1. 

5.1.  Calculations  and  Results 

In  Section  4.1,  the  trial  ground  state 

¥0  =  expUXtrU  } 

P     V  (5-la) 

was  proposed  and  the  variational  parameter  X  was  determined  by  minimizing  the 
ground  state  energy  per  plaquette.  If  I  ^  >  is  a  trial  first  excited  state  of  the  form 
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«F1>=Fl«F0> 


where  F  =  F  -  <F>,  with  F  is  a  linear  combination  of  plaquette  variables,  and 
where  the  expectation  of  F,  <F>,  is  taken  with  respect  to  the  ground  state 
I  *F  >,  such  that  1 4^  >  is  orthogonal  to  I  *F  >, 


<  V  j  I  ¥Q  >  =  0. 


The  variational  parameters  in  F  are  determined  by  minimizing  the  lowest  excited 
state  energy  as  in 


<*Fj  aHl 

aE-  V¥,]  = 


=  {|^<ItrU  >  +  4^(2-trU  )  } 
"S     P  v 


+  (  —  -4^)l^ln<F+F> 
8^        2  d\ 


.  <[E^[E^,F+]]  F>   <F+[  E?,  [  E?,  F  ]  ]  > 

+  -!-!!{  £  +  ^— ^  

164  k  a  <f-p>  <p+p> 


107 


<[E?,  F+][E*,F]> 
-2  *  __k   j 

<F  F>  (5-2) 

The  derivation  of  eq.  (5-2)  is  shown  in  Appendix  B.  The  term  in  the  first  brace 
bracket  is  just  the  variational  ground  state  energy.  Since  the  glueball  mass .  Mg, 
is  defined  as  the  energy  difference  between  the  lowest  excited  state  and  the 
ground  state,  eq.  (5-2)  gives  the  glueball  mass  as 


Mga  =  a(E1-E0) 


=  (— -4£)  j-^-ln  <F+F> 
8^         2  dX 


j  <[  E*  [  E*.  F+ ]  ]  F>   <F  [E^[E^F]]> 

-f-    S  S    {   —  


<[Ej,F+][Ej,F]> 

-2  * — jr1*  } 

<F  F>  (5-3) 


Eq.  (5-3)  can  be  simplified  greatly  by  assuming  the  trial  ground  state, 
eq.  (5- la),  is  exact,  then  the  glueball  mass  is(^l,43) 


l<[F+,[aH,F]]> 

Mga"2  ^  

<F  F> 


108 


<(  E*  F+)  ( E*  F  )> 
=  .  J- IE  — — — £  

k  3         <F+F>  (5-4) 
Here  the  second  equality  is  obtained  by  using  the  integration  by  parts.  The 
results  of  eqs.  (5-3)  and  (5-4)  should  agree  with  each  other  in  the  strong- 
coupling  region  because  the  variational  vacuum  becomes  the  exact  ground 
state  in  that  region.  However,  different  results  should  be  expected  in  the 
intermediate  coupling  region.  While  eq.  (5-4)  is  accepted  commonly  by 
most  theorists(37,40,43)  as  the  starting  point  for  calculation  of  glueball 
mass,  in  this  dissertation  both  expressions,  eqs.  (5-3)  and  (5-4),  will  be 
used  simultaneously  to  perform  the  calculations  of  glueball  mass.  It  is  hoped 
that  by  comparing  the  two  results  some  information  about  the  physical 

system  will  be  obtained. 

To  choose  a  projector  F,  notice  that  ^r—  is  the  correlation  length  in  the 

s 

unit  of  lattice  spacing.  As  determined  by  the  recent  Monte  Carlo  calculations 
of  glueball  mass(37)  for  SU(N)  lattice  gauge  theories  with  N=3, 4,  5,  and 

6,  the  estimated  values  of  -rr—  in  the  intermediate  coupling  region  are 

expected  to  be  less  than  1  for  SU(2)  lattice  gauge  theory.  The  projector,  F, 
therefore,  is  chosen  to  include  up  to  two  adjacent  plaquette  correlations 

F  =  7S={Pl|*l(P)  +  P2S  *2  <P2)  ♦  P3  X  *3  (P3)  ) 

V  P        v  v2  v3  (5.5) 
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where  Pj  are  variational  parameters,  ^(p)  is  the  trace  of  a  plaquette, 
is  the  trace  of  an  extended  plaquette  containing  two  adjacent  plaquettes  lying 
in  the  same  plane,  and  <t>^  (P3)  is  the  trace  of  an  extended  plaquette 
containing  two  adjacent  plaquettes  lying  in  two  different  planes.  These  three 
operators  are  shown  in  Figure  24. 


(c) 

Figure  24.  The  operators  used  in  the  projector  F 

a)  An  elementary  plaquette;  b)  Two  adjacent  plaquettes 
lying  in  the  same  plane;  c)  Two  adjacent  plaquettes  lying 
in  two  different  planes. 


By  using  the  projector  of  eq.  (5-5),  eq.  (5-3)  becomes 


1 


{ 


8         I  P  D.  .  p.  8 
i j=l  1    i.J  J 


p.*N.  .  p.  +-V  .  I    P*Q.N.  .p. 
8  ij=i  Ki    10  KJ    16  i,j=i  *i  ^i,j   10  KJ 


+  (^-4^) 


3 


ij=l 


p*  [ii-D..] 


(5-6) 
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and  eq.  (5-4)  becomes 


.* 


.?  ,  p.  N..  p. 

Mga^=  4  "i — ; — 


M=l    1    y  J  (5-7) 


where 


Ng  =  -(1  +  \2+35i,3><if^<(°»^?^)» 


(5-8) 


D.  .  =  (1+8.  +3  8.  ){  «()*(0)  Z  (|>.(p)>  -  «J)*(0)>  <I  (|>.(p)>  } 
1,J  1,2        1,3       Yi v    p  \p'       Yi v  '     p  Yj 


and 


(5-9) 


15 

15  n 

6 

2 

2 

15 

2 

9 

9 

15 

2 

9 

9  . 

(5-10) 


denotes  the  trace  of  a  itn  type  of  plaquette  at  a  fixed  position  and  the 
overall  factor  (1  +  8.  _  +  38.  J  in  N..  and  D..  results  from  counting  the 

1,Z        l,j         rj  rj 
number  of  the  extended  plaquettes  on  a  lattice  in  the  unit  of  the  total  number 

of  plaquettes  Np.  The  minimization  condition,  _<L  (M  a£)  =  0  >  yields 

3p.  & 


detl  }N  +  1LQD  +  (f-4^2)I|_D.(MgaOD|  =0 


(5-11) 


Ill 


foreq.  (5-6)  and(43) 


detl^N-^M  a)DI  =  0 


(5-12) 


for  eq.  (5-7)  for  the  nontrivial  Pi.  The  solutions  of  eq.  (5-11)  and  eq.  (5-12) 
will  give  the  glueball  mass  for  the  cases  of  eqs.  (5-3)  and  (5-4),  respectively. 


(|>.  (o)  and  <J>j(p)  in  eq.  (5-8),  contributions  to  the  matrix  elements  of  N  then 

come  from  adjacent  and  overlapped  diagrams  formed  by  any  two  plaquettes  of 

various  types  in  eq.(5-5).  Using  the  Gross-Witten  trick,(21)  it  can  be  shown 

that  the  contribution  from  the  adjacent  diagrams  vanishes.  To  see  this,  consider 

an  adjacent  diagram  formed  by  two  plaquettes,  U    and  U    ,  of  various  types 

pl  p2 

in  eq.  (5-5)  with  one  common  link.  Then 


Since  the  color  electric  field  operator,  E^,  acts  only  on  the  common  links  of 


<H(^Pi)(B^,,0P2)>"<J,,(f<,Pi),,(TV> 


(5-13) 


where  the  second  equality  is  followed  by  the  identity 


i.    8. .  -  8. .  8. 
i,m  j,k     i,j  k,m 


) 


(5-14) 
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and  the  result  of  eq.  (5-13)  vanishes,  through  eq.  (4-26)  which  is  obtained  by 

using  the  Gross- Witten  trick.(21)  Therefore,  for  the  matrix  elements  of  N,  only 
the  overlapped  diagrams  of  the  matrix  elements  Ny,  which  are  shown  in 

Tables  5  through  8,  must  be  considered. 

The  matrix  elements  of  D  defined  in  eq.  (5-9)  signify  the  correlation  effects. 

Although  the  variational  vacuum  is  a  product  of  the  single-plaquette 

wavefunction,  correlations  will  be  involved  through  the  lattice  Bianchi  identities 

even  for  disconnected  diagrams,  such  as  two  separated  plaquettes  without  a 

common  link.  Thus  the  calculations  become  complicated.  To  manage  the 

calculations  of  the  matrix  elements  of  D,  the  following  approximations  are  used: 

In  the  zeroth  order,  the  correlation  between  two  extended  plaquettes  (i,j  =  2  and 

3)  is  approximated  by  the  correlation  of  an  extended  and  an  elementary  plaquette. 

That  is,  for  D22  the  approximation  yields 


D22  =  2  {  <  <>2(o)  Z  «>2  (P)  > "  <  Vo)  >  <  S  h  (P)  >  1 


4  < 4^(0)  >  {  < <j>2(o)  I ^(p)  >  -  < <|>2(o)  ><  I  0j(p)  >  } 


=  4<<))1(0)>{}|^<(()2(0)>} 


(0) 

22  (5- 15a) 


where  the  Gross-Witten  trick  was  used  to  get  the  second  equality  and 
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Table  5.     The  overlapped  diagrams  contributed  to  N  jj  for  i  =  1,  2  and  3. 


Matrix  Element 

Type  of  Diagrams 

Number  of  Diagrams 

Nll 

1 

N12 

4 

N13 

A 
//  / 

8 

Table  6.   The  overlapped  diagrams  contributed  to  N  22  • 


Type  of  Diagrams 

Number  of  Diagrams 

2 

4 

8 

Table  7.  The  overlapped  diagrams  contributed  to  N  23  • 


Type  of  Diagrams 

Number  of  Diagrams 

A 

//  / 

4 

/  / 

V 

16 

/ 

/ 

// 

8 
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Table  8.    The  overlapped  diagrams  contributed  to  N33 . 


Type  of  Diagrams 

Number  of  Diagrams 

// 

k 

V 

4 

A  1 

r, 

1 

/ 

8 

// 

A 
I 

16 

// 

* 

/ 
/ 

8 

/ 

V 

A 
/ 

8 

A 
/ 

16 
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the  superscript  on  denotes  the  zerotn  order  approximation.  In  general,  the 
zerom  order  approximation  gives 


Df=I(1+S.2  +  3  8.3)(  1  +  5.2  +  38.3) 

<V0)>,i|-<«.(0)>  +  I|:<*j(0)» 

fori*landj*l  (5-15b) 

To  improve  the  zero**1  order  approximation,  notice  that  the  diagrams  contributed 
to  the  correlations  can  be  divided  into  three  types:  Disconnected  (no  common 
links),  connected,  and  overlapped.  The  dominant  contributions  come  from  the 
overlapped  diagrams.  Therefore,  to  improve  the  approximation  the 
contributions  of  the  overlapped  diagrams  are  calculated  exactly  and  the 
difference  between  the  results  of  the  exact  calculation  and  the  zero**1  order 
approximation  for  the  overlapped  diagrams  is  found.  Finally,  the  difference  is 
added  to  eq.  (5- 15b)  to  give 


D..-.D(0)+D<1) 
y      y       y  for  i*  land  j*l  (5-16) 
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Table  9.  The  matrix  elements  D  „  for  i  *  1  and  i  *  1 . 

■J 


4i2)=(i[<oo>-<nn>2i 


♦s  [<□[□□  >-<nn>2n 
-[*<□  >[<□□  >-<  □  ><□□>]) 


+  4[< 


nj/s/i  ><ljO  ><□□  >i 

>■< ><□□>]> 

(8<  □  >[<  >-<  £7  ><  £jQ  >]) 


D33  =  »< 


>L  ]+2[< 


+  2[< 


>2] 


>2n 


>]} 
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where       are  shown  in  Table  9.  The  matrix  elements  of  D  between 

y 

elementary  plaquettes  and  extended  ones,  i.e.,  Dil  for  i  =  1, 2  and  3,  can  be 
obtained  exactly 


Dil=(1+5i,2  +  35i,3){2-i^i(°)>}- 


(5-17) 


Now  all  the  matrix  elements  of  N  and  D  are  expressed  in  terms  of  the 
expectation  values  of  some  physical  observables  with  respect  to  the  trial  ground 
state  wavefunction,  eq.  (5- la),  and  the  method  developed  in  Chapter  3.2  can  be 
used  to  calculate  these  expectation  values.  Once  the  matrix  elements  of  N  and  D 
are  obtained,  the  glueball  mass  can  be  found  by  solving  eqs.  (5-11)  and  (5-12). 
The  resultant  glueball  mass  are  shown  in  Figure  25.  The  squared  dots  represent 
the  result  of  eq.  (5-3),  the  triangular  dots  are  those  from  eq.  (5-4),  and  the  solid 
lines  are  the  expected  asymptotic  scaling  lines.  It  can  be  seen  from  Figure  25 
that  both  results  exhibit  an  obvious  tendency  toward  scaling,  but  the  glueball 
mass  in  this  calculation  does  not  scale  exactly  along  the  expected  asymptotic 
scaling  line.  By  comparing  the  two  curves  in  Figure  25,  two  remarks  can  be 
made: 

(1)  The  two  results  coincide  with  each  other  for  ^  <  0.1 ,  which  implies 
that  the  variational  vacuum  becomes  the  exact  ground  state  for  £  <  0. 1  . 

(2)  In  the  intermediate  region  there  is  a  significant  difference  between  the 
two  results.  In  particular,  the  curve  falls  down  during  the  range  of 

0.22  <  £  <  0.30  for  the  result  of  eq.  (5-6)  and  0.27  <  £  <  0.35  for  the  result  of 
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eq.  (5-7).  This  implies  that  the  trial  ground  state  wavefunction,  eq.  (5- la),  is 
not  good  enough  to  assume  it  to  be  exact 
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Figure  25.  The  result  of  glueball  mass.  The  squared  dots  represent 
the  result  of  eq.  (5-3),  the  triangular  dots  that  from 
eq.  (5-4),  and  the  solid  lines  are  the  expected  asymptotic 
scaling  lines. 


5.2.  Discussions  and  Comments 

One  intuitively  derived  suspicion  to  explain  why  the  resultant  glueball  mass 
does  not  scale  exacdy  along  the  expected  asymptotic  scaling  line  for  a  certain 
range  of  coupling  constant  is  the  following:  The  correlations  contained  in  the 
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lowest  excited  trial  state  are  not  strong  enough  to  yield  the  necessary  collective 
mode.  The  extended  plaquettes  contained  in  the  projector  are  all  six-link 
plaquettes.  Recent  results  obtained  by  Long,  Robson,  and  Chin(^)  for  SU(3) 
glueball  mass  indicate  that  the  six-link  twisted  plaquette  (see  Figure  26)  may 

A 

i 
i 
i 
i 

J- 

4 


Figure  26.  A  six-link  twisted  plaquette. 

play  a  more  significant  role  than  those  included  in  this  calculation. 

Since  glueball  mass  is  a  quantity  determined  through  the  collective  mode, 
the  scaling  window  is  expected  to  be  very  small.  However,  it  was  pointed  out 
by  Chin  and  Karliner(45)  that,  in  general,  the  ratio 

Mg 
To 

as  a  function  of  coupling  constant  should  behave  like  the  one  depicted  in  Figure 
27:  The  ratio  R  starts  from  the  strong-coupling  region  (no  scaling),  exhibits  a 
scaling  window,  and  eventually  diverges.  Furthermore,  the  scaling  window 
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Figure  27.  The  behaviors  of  R  in  different  coupling  regions. 
[This  figure  is  taken  from  ref.  (45).] 


exhibited  by  the  ratio  R  can  cover  a  wider  range  of  coupling  constants  than  that 
exhibited  by  a  single  physical  quantity  M  or  Jc-        Tne  ratio  R  therefore  is 
computed  and  the  result  is  shown  in  Figure  28.  The  scaling  window  is 
observed  and  the  scaling  behavior  implies  that 


Mg  =  2.867o 

Previous  estimates  of  the  glueball  mass  in  Lagrangian  Monte  Carlo  calculations 
are  M  =(3.5  ±  1.2)  Jo  by  Berg(46)  and  M  =  (  3.1  ±  1.1)  Jo  by  Bhanot 
and  Rebbi.(42)  The  result  obtain  here  agrees  with  those  estimates. 
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Figure  28.  The  ratio  of  glueball  mass  to  the  square  root  of  string  tension. 
The  square  dots  are  obtained  by  using  the  glueball  mass  of 
eq.  (5-3),  and  the  triangular  dots  are  that  from  eq.  (5-4). 


CHAPTER  6 

MONTE  CARLO  SIMULATION  AND  THE  BIASED-SELECTION 

METHOD 

Creutz(4)  successfully  applied  Monte  Carlo  methods  to  the  Lagrangian 
formalism  of  lattice  gauge  theories  to  demonstrate  the  coexistence  of  color 
confinement  and  asymptotic  freedom  in  non-nonabelian  gauge  theories.  Since 
then  the  Monte  Carlo  techniques  have  become  an  indispensable  numerical  tool 
for  investigating  the  nonperturbative  properties  of  lattice  gauge  theories.  The 
Monte  Carlo  methods  are  an  old  technique,  referred  to  as  "  importance 
sampling.  "(10)  For  example,  consider  an  integral  of  the  form 


Jtndu.  looy^cuj 
<o(U)>=    )    k      p  L 

k     k         p  (6-1) 


in  SU(2)  Hamiltonian  lattice  gauge  theory,  where  V(Up)  is  a  known  trial 

ground  state  wavefunction  and  0(Up)  is  a  physical  observable.  If  a  6^  site 
lattice  is  used,  there  are  648  total  linked.  Since  each  link  is  associated  with  a 
group  element  of  SU(2),  there  are  3  independent  variables  per  link,  with  a  total 
of  1944  variables.  Equation  (6-1)  then  becomes  a  1944-dimensional  integral. 
If  a  conventional  numerical  method,  e.g.  M-point  quadrature  formula,  is  used  to 
do  the  integration,  M^944  evaluations  of  the  integrand  are  required.  However 
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in  the  M*944  evaluations  very  few  configurations  contribute  significantly  to 
the  expectation  value  <  0(Up)  >  significantly,  and  much  effort  is  wasted  on  the 
insignificant  configurations.  The  difficulty  in  evaluating 
eq.  (6-1)  is  shifted  to  the  generation  of  important  configurations.  To  state  the 
problem  precisely,  the  following  probability  density 


P(UP)=1  ~2  

p    J  [ndu.iv2  (uj 

k     k         P  (6-2) 


is  defined.  Equation  (6-1)  then  becomes 

<o(up)>  =  j  [  n  duk  ]  0(up)  P(up) 


Now  the  problem  can  be  addressed:  Without  actually  knowing  what  P(Up)  is, 
how  can  an  ensemble  of  configurations  which  are  distributed  according  to 
P(Up)  be  generated?  By  viewing  P(Up)  as  the  stationary  distribution  of  a 
Markov  process,  the  above  problem  can  be  solved  and  the  evaluation  of 
eq.  (6-1)  be  accomplished.  There  are  two  well-known,  widely  used 
computer  algorithms-the  Metropolis  algorithm^'47)  and  the  heat  bath 
algorithm(4>27)  -which  are  adequate  to  process  the  Markov  chain.  These  two 
algorithms  are  discussed  in  Section  6.1. 

The  tremendous  consumption  of  computer  time  is  always  a  problem  in 
performing  Monte  Carlo  calculations,  even  solely  for  the  gluon  sector.  Thus, 
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improvement  of  the  numerical  technique  is  essential  for  a  complete  study  of  the 
theory.  Section  6.2  focuses  on  this  problem  with  yet  another  Monte  Carlo 
technique,  the  Biased-Selection  Monte  Carlo  method. 

(13) 

6.1.  The  Markov  Process 
A  common  problem  in  lattice  gauge  theory  is  the  evaluation  of  an 
integral  of  the  form 


<0(U,  )>  =  J  [  n  dU.  ]  0(U, , )  P  (U, .. ) 

k  k     k        k     s    k  (6  4a) 


where  0(UjP  is  a  physical  observable  and  ^(U^)  is  analogous  to  the 
probability  distribution  function  in  statistical  mechanics.  In  the  Lagrangian 
formalism,  eq.  (6-4a)  is  a  Euclidean  expectation  value  of  00-^)  with 


exp  [  -  S(U_)  ] 

Ps(Uk)  =  -f  E  

k     JtndU,  ]exp[-S(UD)] 

k     k  P  (6-4b) 


where  S(Up)  is  the  lattice  action  defined  in  eq.  (2-5).  In  the  Hamiltonian 
formalism,  eq.  (6-4a)  is  a  variational  estimate  of  a  ground  state  expectation 
value  of  0(U^)  with 


^(U  ) 

k     J[ndU  .UK2  (II) 

k     k  k  (6-4c) 
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where  V(U^)  is  a  known  trial  ground  state  wavefunction.  In  the  Metropolis  and 
the  heat  bath  algorithms,  PS(U, )  *s  regarded  as  the  stationary  distribution  of  a 


probability  W(  U,  ->  U'  )  for  changing  the  link  value  from  U,  to  U'  is  devised 

such  that  an  arbitrary  initial  distribution  PQ  (U^)  will  evolve  into  Ps(UjP 
through  the  Markov  process. 

In  both  algorithms  the  starting  point  is  an  arbitrary  configuration  bearing 
with  an  initial  distribution  PQ  (U^).  A  single-link  variable  then  is  altered  and  a 
new  configuration  {  U,  }  is  generated.  The  new  configuration  is  accepted  or 
rejected  depending  on  whether  the  rules  of  the  algorithm  are  satisfied  or  not. 
After  all  the  link  variables  have  been  sampled,  the  entire  lattice  has  been  "swept 
through"  once  and  a  configuration  {  C.  }  is  established.  Using  the 
configuration  established  previously,  this  process  is  repeated  N  times  to  end  up 
with  N  configurations.  The  expectation  value  <  0(U, )  >  can  then  be  estimated 


Markov  process.  From  knowing  V  (U, )  or  exp  [  -  S(U, )  ],  a  transition 


by 


<0(U  )>  =  -U!  0({C}.)±a 

K         N  1=1  1 


(6-5a) 


where  gz  is  the  variance 


[0({C}.)-<0(Uk)>]2 


(6-5b) 
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However,  there  is  an  underlying  assumption  for  eq.  (6-5b):  The  N 

configurations  are  assumed  to  be  mutually  uncorrelated.  To  avoid  the 

correlation  between  configurations  in  a  Markov  chain,  the  algorithm  can  be 
applied  many  times  on  (C}m  until  a  new,  statistically  independent  member  of 

the  stationary  distribution,  {C}        is  generated.  Since  the  standard  deviation 

o  is  inversely  proportional  to  /N  ,  the  magnitude  of  N  depends  on  how 
accurate  is  the  result  that  is  required.  A  typical  value  of  N  is  about  3000. 
6.1.1.  The  Metropolis  Algorithm 

For  this  algorithm,  one  starts  with  a  configuration  { C } ,  and  a  new 
configuration  is  generated  through  the  following  procedure: 

(47) 

(a)  The  group  element,  U^,  associated  with  the  link,  k,  in  question  is 

replaced  by  Uj^and  all  other  links  are  held  at  their  current  values. 

2 

(b)  The  change  in  the  action  S  (or  in      )  is  computed.  If  AS  is  less  than 

2 

zero  (or  A*F  is  larger  than  zero),  the  change  is  accepted.  If  AS  is  larger  than  zero 
2 

(or  A*P  is  less  than  zero),  the  change  is  accepted  with  the  condition  probability 
v2  (C) 

exp  (-  AS  )  (or  —  ) .  The  transition  probability  w(  U  -»  U'  )  here  is 

assigned  according  to  the  detailed  balance. 

(c)  The  entire  procedure  then  is  repeated  for  the  next  link  until  the 
entire  lattice  is  covered.  This  represents  one  Monte  Carlo  iteration. 

(d)  For  the  next  Monte  Carlo  iteration,  the  previously  established 

configuration  is  used  as  an  initial  configuration  and  then  the  steps  (a)  ~ >  (c)  are 
repeated  to  establish  another  new  configuration. 
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To  show  how  the  procedure  eventually  brings  the  system  into  the 
stationary  probability  distribution  defined  in  eq.  (6-4b)  or  (6-4c), (22,47) 
consider  the  relationship  between  the  ktn  and  the  (k  +  l)tn  sweeps, 

Pk+1(C)  =  S,  Pk  (C)  W(C  -»  C)  +  (1  -  I  W(C  ->  C) )  P  (C) 
{C }  1^1 


=  P  (C)  +  Z  {  P  (C)  W(C->  C)  -  W(C  ->  O  P.  (C)  } 
k     IC)    k  k  (6.6a) 


where  W(C  ->  C)  is  the  probability  for  the  transition  C  ->  C  with 


{Z}W(C-.C)<1 


(6-6b) 


For  a  stationary  distribution, 


P^1(C)  =  P^(C) 


(6-7) 


From  eq.  (6-6a),  this  implies  the  detailed  balance 

P*  (C)  W(C  ->  C)  =  Pf  (C)  W(C  -»  o 
k  ^  (6-8) 

Now  the  transition  probability  in  the  Metropolis  algorithm  is 

WCC-XW1  ifS(C)<S(C) 

exp  [  -  (S(C')  -  S(C»  ]      if  S(C)  >  S(C)  (6.9a) 
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or 


W(C^C)^V(0 


V2(Q  (6-9b) 


if\|^(C)>\|/2(C) 
if\|T2(C)<\|/2(C) 


So  then 

W(C  -*  C) 


=  exp[-(S(C)-S(C))] 


W(C-»  Q      r  (6.10a) 


or 


W(C  ->  C)  _  V2  (C) 
W(C->Q  ^(q 


or 


(6-  10b) 


From  eq.  (6-8),  this  imphes 

P*  (C)  =  N  exp[  -  S(C)  ] 
K  (6-1  la) 


P*(C)  =  Ny2(C) 


(6-1  lb) 


where  N  is  a  constant  which  is  determined  by  the  normalization  condition 


L  P*  (C)  =  1 

lc)  k  (6-llc) 


Therefore,  the  algorithm  provides  the  stationary  probability  distribution  defined 
in  eq.  (6-4b)  or  (6-4c). 

The  Metropolis  algorithm  is  very  easy  to  implement  because  the  explicit 
structure  of  the  group  measure  can  be  ignored.  However,  in  a  weak  coupling 
region,  the  correlation  length  of  the  system  becomes  large  and  most  of  the 
changes  therefore  would  be  routinely  rejected.  The  configurations  thus 
obtained  are  highly  correlated.  This  is  a  practical  problem  in  applying  the 
Metropolis  algorithm. 
6.1.2.  The  Heat  Bath  Algorithm 

Given  the  heat  bath  algorithm  (4,30,48)  ^  starting  with  a  given 
configuration,  a  new  configuration  is  generated  through  the  following 
procedure: 

(a)  Starting  from  a  link  k,  all  of  the  links  which  interact  with  k  have  fixed 
U  matrixes  and  provide  a  thermal  background.  There  are  six  plaquettes 
interacting  with  k  in  the  Lagrangian  formalism,  and  four  plaquettes  in  the 

Hamiltonian  formalism. 

(b)  The  group  element,  U, ,  of  link  k  is  replaced  with  a  new  element,  U'  • 

k  k 

chosen  randomly  from  the  entire  group  with  its  probability  proportional  to  the 
Boltzmann  factor 
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dP(U')  ~  exp[  -S(U')  ]  dU1  (6.12a) 

in  the  Lagrangian  formalism  or  to 

dP(U)  ~  ¥2  (IT)  dU1  (6.12b) 

in  the  Hamiltonian  formalism  S(U')  (or  y2  (U'))  is  evaluated  with  the  link  k 

having  the  value  U'  and  all  other  interacting  links  fixed  at  their  previous 
values. 

(c)  This  procedure  is  applied  to  every  link  in  the  lattice  to  complete  one 
iteration. 

(d)  For  the  next  Monte  Carlo  iteration,  the  previously  established 
configuration  is  taken  as  the  initial  configuration  and  then  the  steps  (a)  to  (c)  are 
followed  to  finish  another  Monte  Carlo  iteration. 

The  heat-bath  algoritnm  also  satisfies  the  detailed  balance.  It  is  more 
effective  than  the  Metropolis  algorithm  because  for  each  iteration  the  system  is 
brought  closer  to  the  stationary  probability  distribution.  However,  the  heat 
bath  is  harder  to  implement  because  it  needs  the  group  measure  explicitly  in 
eq.  (6-12a)  or  (6- 12b).  For  complicated  gauge  groups  like  SU(3),  it  is  tedious 
to  work  out  eq.  (6- 12a)  or  (6-12b).(48)  Thus,  most  Monte  Carlo  simulations  in 
practice  have  turned  to  the  Metropolis  algorithm. 
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6.2.  The  Biased-Selection  Monte  Carlo  Method 
The  Biased-Selection  Monte  Carlo  method  is  a  multi-dimensional  extension 
of  information  sampling.  It  was  introduced  by  Coldwell, 

(13) 

who  used  a 

simple  form  of  this  method  to  calculate  the  partition  function  of  a  hard  disk 
system  in  statistical  mechanics  and  more  recently  applied  this  method  to 
calculate  the  potential  between  two  helium  atoms.O^) 
In  the  Biased-Selection  method,  an  integral  such  as 

I(f)  =  Jdx1dx2-  dxn  f(Xj,  x2,  •••  ,  xn  ) 


is  estimated  by 


Kf)  =  ^.£   1   2W  ±a, 

N  1=1  W.  1 

i  (6- 14a) 


where  each  configuration  { xp  x2,  •  •  • ,  xn  } .  is  a  set  of  values  for 

xp  x2,  •  •  • ,  xn  which  was  selected  with  probability  Wi  relative  to  the  uniform 
and  random  distribution,  and  the  standard  deviation  Q"j  is  given  by  the  usual 

expression 


1 


(6-14b) 
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Expectation  values  such  as 


^      JdXjdXj  -  -  dxn       (xr  x2,       xn)  f(xr  x2,  xn) 
jdx^Xj---  dxn  \\f2  (xr  x2,  xn) 


(6-15) 


are  estimated  by 


.1         {xr  x2,  -,  xn}.  )  f  (  {xr  x2,      xn}.  )  /  W. 

 N  °f 

|  v^({Xl,  x2,-,  x  }  )/W. 

1=1  12        n  1         1  (6-16a) 


with  standard  deviation 

I  [^({x}  )(f({x}.)-<f>)/W.]2 

d=m  -i  !  l_ 


[i  vz({x}  )/w  r 

1=1  1       1  (6- 16b) 


The  standard  deviation  for  <  f  >  is  smaller  than  that  for  a  ratio  of  two  integrals 
since  the  same  configurations  are  used  for  both  the  numerator  and  the 

denominator.  W.  is  determined  by  the  guiding  function  used  to  choose  the 
configurations.  When  the  method  for  selecting  configurations  is  a  Markov 
chain,  y  is  used  as  the  guiding  function  and  eq.  (6- 16a)  reduces  to  an  average 
over  equally  weighted  values  of  f.  Therefore,  the  Markov  chain  is  a  subset  of 
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the  Biased-Selection  method.  But  since  correlations  between  configurations  in 
the  naive  Markov  chain  are  very  high,  eq.  (6- 16b)  is  not  true. 

The  following  section  discusses  in  detail  the  application  of  this  method  to 
the  particular  case  of  calculating  the  ground  state  energy  per  plaquette  for  SU(2) 
Hamiltonian  lattice  gauge  theory.  The  results  demonstrate  not  only  the 
accuracy,  but  also  the  speed,  of  the  method. 
6.2.1.  Algorithms  and  Examples 

As  discussed  in  Section  4.1,  if  the  wavefunction 


is  used  as  the  trial  ground  state  where  X  is  a  variational  parameter,  then  the 
ground  state  energy  per  plaquette  is 


aE0=(^-4$)P(X)  +  85 


(6- 18a) 


where  £  =  — 5— ,  and  P(K)  is  the  expectation  value  of  plaquette 


j[IIdU,  ](LtrU_)v5 


(6- 18b) 
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The  function  P(X)  in  this  form  is  the  function  to  be  evaluated  using  the  Biased- 
Selection  method. 

There  are  two  important  ingredients  in  the  following  algorithm.  The  first  is 
the  guiding  function  <t>  for  the  biased  selection  of  the  configurations.  It  is  not 
necessary  for  the  guiding  function  to  be  the  same  as  Yo  •  1°  most  cases,  a 
single  guiding  function  is  used  to  cover  certain  ranges  of  X  in  ¥o  if  the  physical 

properties  implied  by  Vo  do  not  change  dramatically  over  those  ranges  of  ~k. 

Another  important  nature  of  the  algorithm  ties  in  the  partitions  of  lattice. 
When  a  function  f(x,y)  is  given  approximately  by  f1(x)f2(y),  a  product  of 
probabilities  (which  corresponds  to  the  case  of  using  f(x,y)  to  choose  x  and  y) 
can  be  turned  into  a  sum  of  probabilities  by  using  f  t(x)  to  choose  x  and  f2(y)  to 
choose  y.  In  the  case  of  the  integration  in  eq.  (6- 18b)  where  the  trial 

wavefunction,  eq.  (6-17),  is  a  product  of  single-plaquette  wavefunction,  the 
integration  can  not  be  reduced  to  a  single-plaquette  integration.  As  discussed  in 
Chapter  3.2,  the  correlations  between  various  plaquettes  will  appear  through  the 
effect  of  Bianchi  identities  because  the  Haar  measure  is  defined  in  the  link 
space.  However,  the  results  of  Chapter  3.2  showed  that  the  effects  of  the 
Bianchi  identities  are  very  localized  from  the  strong-coupling  region  down  to 
the  intermediate  coupling  region.  The  wavefunction,  eq.  (6-17),  therefore  can 
be  sampled  in  terms  of  groups  of  plaquettes.  The  larger  the  effects  of  Bianchi 
identities,  the  bigger  will  be  the  number  of  plaquettes  in  a  group.  If  a  group  of 
plaquettes  can  be  sampled  adequately,  then  a  large  step  will  have  been  taken 
toward  sampling  the  full  trial  ground  state  wavefunction. 
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In  the  case  discussed  here,  the  number  of  plaquettes  in  a  group  is  given  as 
one.  Starting  with  a  particular  plaquette,  P,  and  a  guiding  function,  <}>,  a 
configuration  is  generated  by  the  following  procedure: 

(a)  The  group  elements  of  SU(2)  correspond  to  the  unit  vectors  in  four 
dimensions.  The  four  links  of  the  plaquette  P  are  assigned  values  by  randomly 
choosing  values  for  the  four-unit  vectors  associated  with  each  link,  L  times. 

(b)  We  then  form 


with  j  <  L,  and  where  <{>  =  n  (t>p(Up).  Next  a  random  number,  R,  is  chosen 
between  0  and  S(L)  to  pick  up  a  set  of  link  values,  U  (k),  ne  p,  such  that 


(6-19) 


(6-20a) 


The  probability  for  R  being  in  the  interval 


S(k-l)<R<S(k) 


(6-20b) 


is  exactly 


P[  Up(k)]  = 


yupoo) 

{^S(L)} 


(6-21) 
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Therefore,  if  Up(k)  is  chosen  for  plaquette  P,  the  integrand  must  be  divided  by 
P[  Up(k)  ]  to  correct  for  this  biasing. 

(c)  For  the  next  plaquette,  P,  which  has  a  common  link,  n,  with  plaquette 
P,  we  go  through  the  same  procedure  (a)  — >  (b)  except  that  the  common  link  n 
is  held  fixed  in  the  previous  value  obtained  from  sampling  the  plaquette  P. 

(d)  This  procedure  is  repeated  until  the  whole  lattice  is  sampled,  and  a 
Monte  Carlo  iteration  is  established.  For  Np  plaquettes  in  the  whole  lattice,  this 
process  requires  a  total  of  LNp  plaquette  evaluations  and  the  probability  of  this 
configuration  being  selected  is  thus 


np 

n  <t> 
p;=i  pi 

w.= 


1    np  1 
n  ff  s  (L)} 

p.=l   L  Pj 


(6-22) 


The  probability  in  eq.  (6-22)  indicates  that  the  LNp  choices  with  the  partitions 

NF 

of  lattice  allow  concentration  on  the  region  found  only  by  L  random  choices 
without  the  partitions  of  lattice.  Therefore,  a  product  of  probabilities  has  been 
converted  into  a  sum. 

(e)  The  process  (a)  — >  (d)  can  be  repeated  N  times  to  accomplish  N 
iterations.  The  N  configurations  thus  established  are  statistically  independent. 

If  Spj(L)  in  eq.  (6-22)  was  the  same  for  each  configuration,  then  eq.  (6- 16a) 
would  reduce  to  a  sum  over  equally  weighted  values  of  the  integral  just  as  in  the 


case  of  the  Markov  process.  In  addition,  since  each  configuration  (in  the  Biased- 
Selection  method)  is  independent  statistically,  eq.  (6- 16b)  is  valid  whenever  there 
are  enough  configurations  to  estimate  errors  through  such  fluctuations.  This  leads 
to  the  definition  of  the  effective  number 


N  „  =  ■  1-1 


eff       N    a  2 
Z^(i)/\\f 
1=1  u  1 


(6-23) 


which  measures  the  number  of  configurations  contributing  non-zero  amounts  to 
eq.  (6-  16a).  Neff  needs  to  be  on  the  order  of  30  or  more  for  eq.  (6- 16b)  to  be 
valid.  Some  examples  can  now  be  discussed: 

In  2  +  1  dimensions  with  free  boundary  condition,  it  was  shown  in 
Chapter  3.2  of  this  thesis,  that  because  the  Jacobian  of  transforming  from  the 
link  space  to  plaquette  space  is  the  identity,  the  integration  of 
eq.  (6- 18b)  can  be  worked  out  analytically.  The  result  is 


P(X)  =  -i-+2  u 


X  1,(4*.) 

1  (6-24) 


where  ln(4X)  is  the  modified  Bessel  function  of  rank  n  defined  in  eq.  (3-3  lb). 

To  test  the  above  algorithm,  the  Biased-Selection  method  is  used  to  calculate 
p(^)  defined  in  eq.  (6- 18b),  and  the  results  are  compared  with  that  of 
eq.  (6-24).  The  guiding  function  (j>  used  for  the  biased  selection  is 
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<KP)=exp[2pXtrtJp] 


X. 

where  P  =  — — -  and  p  will  be  referred  to  as  the  weighting  parameter.  The 
considered  range  of  the  variational  parameter  0  <  A  <  0.6  has  been  divided  into 
three  equal  subintervals,  and  the  midpoint  of  each  of  these  subintervals  has  been 

chosen  as  the  corresponding  weighting  parameter  in  the  guiding  function.  The 

2 

calculation  was  performed  on  a  6  lattice,  L  was  set  at  20,  and  the  number  of 
the  established  configurations  was  1000.  The  result  of  this  calculation  is  shown 
in  Figure  29  (the  squares)  along  with  the  analytical  solution  (the  solid  line). 
There  are  two  striking  features  of  this  calculation:  The  first  is  that  the  Biased- 
Selection  method  produced  the  whole  analytical  result  for  the  range  of  X 
considered  (the  error  is  within  the  size  of  the  squares  in  the  figure).  The  second 
striking  feature  of  this  calculation  is  that  the  required  CPU  time  was  20  minutes 

on  a  VAX  1 1-750.  The  need  to  divide  the  range  of  X  into  subintervals  arises 
from  the  fact  that  if  the  whole  range  of  X.  were  used  for  the  calculation  with  p  = 

0.3  (i.e.,  one  guiding  function),  it  would  result  in  concentrating  the  sampling  at 
higher  values  of  the  integrand  near  X  =  0  (due  to  the  wider  spreading  of  Vo  ^) 
and  it  would  underestimate  the  integral  near  X  =  0.6  (due  to  the  narrower 
spreading  of  Vo  (>-)). 

3 

In  3  +  1  dimensions,  the  calculation  was  performed  on  a  4  lattice  by  using 
the  guiding  function  of  eq.  (6-25).  In  this  calculation,  the  considered  range  of 
the  variation  parameter,  0  <  X  <  0.6,  was  divided  into  five  equally  spaced 


intervals.  L  was  set  at  200  and  the  number  of  the  established  configurations 
was  3000.  This  calculation  took  30  minutes  of  CPU  time  on  an  IBM  3090,  and 
the  result  is  shown  in  Figure  30.  The  error  bars  are  too  small  to  draw  in  the 
figure.  The  result  obtained  here  is  in  best  agreement  with  the  Monte  Carlo 
result  obtained  by  Chin  et  al. 

(38) 

6.2.2.  Discussion  and  Comments 

By  comparing  the  two  calculations  presented  above,  one  can  see  that  the 
CPU  time  required  in  the  case  of  3  +  1  dimensions  is  considerably  larger  than 
that  in  2  +  1  dimensions.  The  reason  can  be  traced  back  to  the  extra  degrees  of 
freedom  in  3  +  1  dimensions.  As  pointed  out  in  Chapter  2.2,  the  correlation 
effect  of  a  lattice  system  depends  on  the  dimensionality  of  the  lattice.  Higher 
dimensions  imply  more  correlations.  It  is  expected,  therefore,  that  the 
correlation  effect  will  be  more  profound  in  3  +  1  dimensions  than  that  in  2  +  1 
dimensions.  The  picture  is  more  transparent  if  the  calculations  defined  in  the 
link  space  are  transformed  to  that  defined  in  the  plaquette  space.  In  the 
Hamiltonian  formalism,  the  Jacobian  for  this  transformation  is  the  identity  in  2 
+  1  dimensions  and  the  Bianchi  identities  in  3  +  1  dimensions.  Therefore,  the 
links  in  the  whole  lattice  are  correlated  in  3  +  1  dimensions.  However,  when 
partitioning  the  lattice  the  number  of  plaquettes  in  a  group  was  chosen  to  be  one 
in  the  above  calculations.  This  results  in  an  inability  to  perform  the  biased- 
selection  for  33%  of  the  total  number  of  plaquettes,  those  which  are  fixed 
during  the  procedure  for  the  biased  selection  of  the  rest  of  the  plaquettes. 

To  improve  the  above  situation,  the  intuitively  derived  step  would  be  to 
increase  the  number  of  plaquettes  in  a  group  which  is  the  basic  unit  for  the 
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partitions  of  the  lattice.  The  calculation,  while  including  three  plaquettes  in  a 

group  for  the  biased  selection  has  been  performed.  The  result  indicates  not  only 
that  the  N  ^  increases  40%  more  than  in  the  previous  case,  but  the  interval  of  X 

for  the  same  weighting  factor  also  increases. 

The  Biased-Selection  method  is  more  complex  than  the  Markov  chain.  For 
quantum  mechanical  calculations,  the  ability  to  remove  the  Coulomb 

singularities  (in  atomic  calculation)  in  E  ( I  xr  x2,  •  •  • ,  xn  I )  with  a  guiding 

function  other  than  Xj/2  ,  along  with  the  need  for  an  accurate     in  local  energy 

calculation,  provides  a  very  good  motive  for  tolerating  the  extra  complexity. 

Finally,  two  additional  remarks  can  be  made  about  the  applications  of  the 
Biased-Selection  method: 

(a)  The  same  configurations  are  used  to  find  the  expectation  values  for  the 
different  values  of  the  variational  parameter  in  the  Biased-Selection  method. 
Therefore,  the  errors  in  the  difference  between  two  expectation  values  are  about 
on  the  same  order  of  magnitude  as  the  errors  in  the  expectation  values.  The 
error  in  the  derivative  of  a  quantity  is  consequently  reduced  greatly  compared 
with  the  case  of  the  Markov  chain. 


Figure  29.    The  expectation  value  of  a  plaquette  in  2  +  1 
dimensions.  The  solid  line  is  the  analytical 
result,  and  the  squares  are  the  Monte  Carlo  result 
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Figure  30.    The  plaquette  expectation  value  in  3  +  1  dimensions. 
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(b)  The  Biased-Selection  method  is  particularly  suitable  for  the  calculations 
of  QCD  at  finite  temperature  and  for  a  full  study  of  the  phase  transitions  that 
occur.  Configurations  can  (without  much  difficulty)  be  generated  on  both  sides 
of  the  phase  transition,  and  included  in  the  sums  in  eq.  (6- 16a).  For  a  given  set 
of  parameters,  an  expectation  value 

t  W?f/W.  +iw?f/W. 
^  _  i=i  i       iu  j=iTi  ic 

I  V?/W.  +£  w?/W. 

i=i  1      lu  i=i  1      1C  (6-26) 

where  subscripts  c  and  u  refer  to  confined  and  unconfined  phases  respectively, 

automatically  selects  the  important  configurations  from  the  appropriate  region. 
Since  in  the  Markov  chain  the  W.  's  are  not  calculated,  these  comparisons  are 
much  more  difficult. 


CHAPTER  7 
CONCLUSION 


Lattice  group  theories  have  offered  a  powerful  approach  for  studies  of 
quark  confinement,  hadronic  spectroscopy,  and  ultrarelativistic  heavy-ion 
collisions,  which  are  closely  related  to  the  nonperturbative  properties  of  the 
-underlying  continuum  theory.  In  this  dissertation  an  attemp  has  been  made  to 
study  extensively  the  computation  methods  for  lattice  gauge  theories  and  the 
physical  properties  of  SU(2)  Hamiltonian  lattice  gauge  theory. 

In  Chapter  2,  the  formulation  of  SU(N)  lattice  gauge  theory  for  both  the 
Lagrangian  and  the  Hamiltonian  formalisms  is  established.  Some  conceptual 
issues  regarding  the  formulations  of  lattice  gauge  theories  are  discussed  in 
detail.  A  way  to  extract  the  continuum  physics  from  a  lattice  calculation  also  is 
set  up.  Chapter  3  is  devoted  to  discussions  of  the  analytical  calculation 
methods.  The  difficulty  caused  by  the  non-physical  singularities  in  the 
calculations  of  strong-coupling  expansion  is  discussed.  It  was  pointed  out  that 
the  variational  calculations  can  circumvent  the  difficulty.  An  analytical  method 
then  was  developed  to  accomplish  the  variational  calculations  of  the  static  quark 
potential  and  the  glueball  mass  for  SU(2)  Hamiltonian  lattice  gauge  theory.  By 
examining  the  analytical  method  in  detail,  the  nature  and  general  rules  for 
applications  of  the  method  were  established.  In  Chapter  4,  the  analytical 
method  was  applied  to  the  variational  calculation  of  static  quark  potential.  From 
the  static  quark  potential,  it  was  possible  to  locate  the  rough  transition  point  and 


144 


145 


extract  the  string  tension  which  is  consistent  with  the  expected  asymptotic 
scaling  line.  The  glueball  mass  for  SU(2)  Hamiltonian  lattice  gauge  theory 
was  calculated  in  Chapter  5  and  the  glueball  mass  in  continuum  physics 
was  obtained  from  the  scaling  window  of  the  ratio  of  glueball  mass  to  the 
square  root  of  string  tension.  In  chapter  6,  a  new  computer  algorithm—the 
Biased-Selection  method— was  developed  to  perform  Monte  Carlo  simulations. 
The  numerical  simulations  performed  in  this  dissertation  by  using  the 
Biased-Selection  method  demonstrated  that  the  new  algorithm  is  accurate 
and  also  is  reasonable  in  its  consumption  of  computer  time. 

There  are  two  main  contributions  contained  in  this  dissertation:  First,  the 
establishment  of  new  computation  methods  and,  second,  the  results  obtained  in 
the  variational  study  of  SU(2)  Hamiltonian  lattice  gauge  theory.  From  a 
methodological  viewpoint,  two  new  computation  methods  are  established 
through  detailed  examinations  and  practical  demonstrations.  One  is  the 
analytical  method  which  permits  performance  of  the  variational  calculations 
analytically.  The  other  is  a  new  computer  algorithm,  the  Biased-Selection 
method,  which  enables  performance  of  the  Monte  Carlo  simulations  correctly 
and  efficiently. 

The  results  achieved  in  the  variational  study  of  SU(2)  Hamiltonian  lattice 
gauge  theory  are  the  following: 

(a)  The  roughening  transition  point  is  located  at  £c  =  0.26. 

(b)  The  scaling  string  tension  is  obtained,  and  thus  the  scale  parameter 


is  determined 
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A„  =  0.009  To 

rl 

(c)  The  scaling  window  for  the  ratio  of  glueball  mass  to  the  square  root  of 
string  tension  is  observed,  and  it  yields 

M  =  2.86  Jo 

These  results  are  helpful  in  understanding  the  underlying  physics  of  SU(2) 
Hamiltonian  lattice  gauge  theory  and  the  physics  pictures  thus  established  can 
possibly  extend  to  the  case  of  SU(3). 

At  present,  the  study  of  lattice  gauge  theory  via  the  Hamiltonian  technique 
is  still  under  development  and  various  problems  of  evaluating  physical 
observables  largely  remain  open.  In  particular,  the  fermion  problem  is  a  very 
challenging  one.  Fermions  are  not  included  in  any  calculation  performed  in  this 
dissertation.  This  is  mainly  due  to  some  surprisingly  difficult  problems  met  in 
the  formulation  of  lattice  fermions  that  have  resisted  resolution.  It  is  generally 
believed  that  some  conceptual  breakthrough  is  needed  to  obtain  a  consistent 
formulation  of  lattice  fermions.  As  mentioned  in  Section  2. 1.2  ,  troubles 
encountered  in  trying  to  consistently  include  fermions  in  the  formulation  can  be 
traced  back  to  the  distorted  space-time  symmetry  of  the  Poincare  group  in  the 
lattice  formulation.  To  pin  down  the  problem,  the  first  step  is  to  formulate 
fermions  in  the  random  lattice  where  the  space-time  structure  is  more  closed  to 
the  reality  than  that  of  the  hypercubic  lattice.  Additional  research  activity 
currently  being  pursued  by  the  author  is  concentrated  on  this  direction. 
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APPENDIX  A 
THE  DERIVATION  OF  EQ.  (4-14) 


If 


li>^+(nx,ny)F.¥(0,  0)lVo> 

and 


<UplVo>  =  exp{|-} 


with  S  =  2X  Z  trUp  ,  then  the  matrix  element  of  between  states  I  i  >  and  I  j  >  for 
the  color-electric  term  in  the  Hamiltonian  is 


< i  I  Z Z Ef Ef  lj> 
k  a   k  k  J 


=  J  [Fldiy  {  exp(|)[E*E*]exp(|-)  } 

=  -J[ndUk]{(E^exp(|))(E^exp(|))F:F. 


+  exp(|)(E*FjH)F.(E£exp(|-)) 


+(E^exp(|-))F:(E^Fi)exp(|) 
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+  exp(f)(E^F;)(EakF.)exp(f)) 


where  the  second  equality  is  obtained  by  using  the  integration  by  parts,  and  the 
summation  in  the  indices,  k  and  a,  is  implied.  The  first  term  in  the  bracket  of 
eq.(A-l)  can  be  rewritten  as 


-  J  [  n  dUk  ]  {  (E*  exp(  | ) )  (E*  exp(  |- ) )  F+  F.  } 

=  -  <VQI  ( E*  lnVoX  E*  lnVo)  F+F.  I  W 

=  l<Vol(E^E^ln¥o)F:F.IV(T> 

+  lJ[jldUk]exp(f)(EjFr)Fi(Ejexp(f)) 


+ 1 J  [  n  dU  ]  exp(  A )  F+  (E?  F  )  (E?  exp(  £ ) ) 
z      k     *         2JK1K        2  (A-2) 


Substituting  eq.  (A-2)  into  eq.  (A-l),  it  yields 

<VoIF/EkEkFi'%> 

=  l<¥ol{[E^[E^,S]]}FjhF.I¥o> 
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-litndUk]exp<f)(E^F;)F.(Eakexp<f)) 

-ij[pdUk]exp<f)FjXF.)(Ekexp<f)) 
-J[5dUk]exp(f)(EkFr)(^F.)«p(f)  (a, 

-Consider  the  following  term  in  eq.  (A-3) 

J  [  n  dUk  ]  exp(  | )  (E*  Fp  F.  (E*  exp(  |- ) ) 

=  -  J  [  n  dUk  ]  (E*  exp(  | ) )  (E*  Fp  F.  exp(  | ) 

-  J  [  n  dUk  ]  exp(  | )  (E*  E*  Fp  F.  exp(  |- ) 

-  J  [  IT  dUk  ]  exp(  | )  (Ek  Fp  (Ek  F.)  exp(  |- ) 

i.  e.,  by  using  the  integration  by  parts,  one  can  show  that 
J  [  n  dUk  ]  exp(  | )  (E*  Fp  F.  (Ek  exp(  | ) ) 

=  -|J  [ndUk]exp(|-)(EkEk  Fp  F.  exp(  | ) 
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-iJ[ndUk]exp(f)(EakF;)(E^.)exp(f) 
Using  a  similar  procedure  of  obtaining  eq.  (A-4),  it  yields 

J  [  n  dUk  ]  exp(  | )  Ft  (Ek  F. )  (Ek  exp(  | )) 


iJ[ndUk]exp(f)F:(EkEkF.)exp(f) 


-Ij[ndUk]exp(f)(EkFr)(EkF.)exp(f) 

Substituting  eqs.  (A-4)  and  (A-5)  into  eq.  (A-3),  the  result  is 

<X|/oIFJEkEkFilx»/o> 

=  i-<Y0l{[E*,[E*,  S]]}F:F.IVo> 


+  l<VoIFj[Ej,[Ej,Fi]]lV(f> 
+  i<Vol[E^[E^Fjh]]F.IVo> 
-i^ltE^FptE^F.]  \%> 


(A-4) 


(A-5) 


(A-6) 


Equation.  (4-14)  can  be  obtained  by  working  out  each  term  in  eq.(A-6)  explicitly. 


APPENDIX  B 
THE  DEPJVATION  OF  EQ.  (5-2) 


If 


I  v,>sFI  vn>  =  (F  -  <F>)  I  wn> 

where  F  is  a  function  of  plaquette  variables,  and  the  expectation  value  of  F, 
<F>,  is  taken  with  respect  to  I  y0  >,  then 

_  ^laHl  yt>  _  ■qfjIaHl  y^/o^l 
i"    OKjIv^     "    oir^V^/^o'V^  (B-2) 

and 

<¥i'Vi> 

 \    1    =<F  F> 

<Vo'V0>  (B.3) 

with 


aH  =  — £  I  EfEf  +  44  2  ( 2  -  trlL  ) 
4^kakk     ^  p  v         P ' 


The  contribution  of  color-electric  term  to  eq.  (B-2)  is 


(B-4) 
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a^a 


<¥,  1 1 1 E,  E.  I  >|f,> 
1     Yi  k  a   k  k  Yi 

45"  c^iy^ 


^  {  <E  trU„>  +  4-^-  In  <F+F>  } 

8$    P    p    2  ax. 


+  -JL— i—  Z  J  {  <[  Ej,  [  E*  F+]  F>  +  <F  [  E*  [  E*  F  ]> 
165  <F  F>  k 


-  2  <  [  E? ,  F+  ]  [  Ef ,  F  ]  >  } 

k      1  L  k'    1    '  (B-5) 


by  using  eq.  (A-6)  in  Appendix  A  directly.  It  is  straightforward  to  calculate  the 
contribution  of  color-magnetic  term  to  eq.  (B-2),  and  the  result  is 


4^ 


oiTj  1 1  ( 2  -  trUp )  I  Vl> 


=  45(<Z(2-trU  )>-l^-ln<F+F>) 

P  P       2  BX  (B-6) 


Combining  eqs.  (B-5)  and  (B-6),  the  result  of  eq.  (5-2)  is  then  obtained. 
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